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The evolution of the Universe is the ultimate laboratory to study fundamental physics across 
energy scales that span about 25 orders of magnitude: from the grand unification scale through 
particle and nuclear physics scales down to the scale of atomic physics. The standard models of 
cosmology and particle physics provide the basic understanding of the early and present Universe 
and predict a series of phase transitions that occurred in succession during the expansion and 
cooling history of the Universe. We survey these phase transitions, highlighting the equilibrium and 
non-equilibrium effects as well as their observational and cosmological consequences. We discuss the 
current theoretical and experimental programs to study phase transitions in QCD and nuclear matter 
in accelerators along with the new results on novel states of matter as well as on multifragmentation 
in nuclear matter. A critical assessment of similarities and differences between the conditions in 
the early universe and those in ultrarelativistic heavy ion collisions is presented. Cosmological 
observations and accelerator experiments are converging towards an unprecedented understanding 
of the early and present Universe. 
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I. INTRODUCTION 



The current knowledge of the early and present Universe is summarized by the standard models of cosmology 
and of particle physics. Their symbiosis provides an unprecedented understanding of the evolution of the Universe 
solidly based on a wealth of observations and experiments. Particle and nuclear physics in a field theory context 
provide the fundamental building blocks which when combined with general relativity and statistical mechanics yield 
a description from the early inflationary stage through a detailed thermal history of the Universe, to the formation 
of large scale structure, galaxies and stars. During the last decade a large body of observational data has provided 
a strong evidence in support of theoretical ideas of an early stage of inflation, during which the visible size of the 
Universe grows exponentially. After this brief, but explosive period of inflation followed by deccelerated expansion 
and cooling, the Universe succesively visits the different energy scales at which particle and nuclear physics predict 
symmetry breaking phase transitions. The goals of this article are to review the equilibrium and non-equilibrium 
aspects of these phase transitions, their cosmological imprints and the current theoretical and experimental efforts to 
study them with accelerator experiments. The article begins with an account of both standard models, setting the 
stage for a discussion of phase transitions in particle and nuclear physics. A brief self-contained excursion of cosmology 
beginning from inflation follows the early history of the Universe, visiting the time-marks at which particle physics 
predicts phase transitions and exploring their potential consequences. The last phase transition(s) of the standard 
model of particle physics took place when the Universe was about 1G~^ sees, old, they are the deconfinement- 
confinement and chiral phase transitions predicted by Quantum Chromodynamics (QCD). The article presents a 
summary of the theoretical efforts and the experimental programs in several accelerator facilities to study the QCD 
phase transitions as well as static and dynamical aspects of phase transitions in cold nuclear matter. 



II. THE STANDARD MODELS: 



The modern understanding of the early and present Universe hinges upon two standard models: the standard 
model of cosmology and the standard model of particle physics^]- 10]. Both have passed stringent observational and 
experimental tests. A wealth of cosmological data is providing confirmation of theoretical ideas in early Universe 
cosmology. Measurements of the temperature anisotropy of the Cosmic Microwave Background radiation (CMB) by 
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satellite, balloon borne and earth based observations, large scale structure surveys, Lyman a forest, cluster abundance, 
weak lensing and measurements of the cosmological expansion and its acceleration by Type la supernovae searches, 
combined with high precision measurements of light element abundances provide an impressive body of high quality 
data that yield an unprecedented understanding of cosmology. 



A. Observational ingredients 

The main observational pillars that support the standard model of cosmology are[ll|-[lol|: 

• Homogeneity and isotropy: on scales larger than > 100 Mpc the Universe looks homogeneous and isotropic. 
This is confirmed by large scale surveys and by the almost isotropy of the CMB. 

• The Hubble expansion: objects that are separated by a comoving distance d recede from each other with a 
velocity v = H d, with H the Hubble parameter, whose value today is Hq « 72 km/s/Mpc. The Hubble 
expansion law determines the size of our causal horizon: objects separated by a comoving distance 

dff = 3000/1^1 Mpc , /i = i?o/100km/s/Mpc (1) 

recede from each other at the speed of light and are therefore causally disconnected. 

• The Cosmic Microwave Background radiation (CMB): a bath of thermal photons with an almost perfect Planck 
distribution at a temperature Tn = 2.725±0.001 K. Temperature anisotropics AT/Tq ~ 10^^ were first measured 
in 1992 by the COBE satellite [ll| , their detection represents a triumph for cosmology. This small temperature 
anisotropy, whose existence is predicted by cosmological models, provide the clue to the origin of structure. It 
is an important confirmation of theories of the early Universe. 

• The abundance of light elements: observations of the abundance of elements in low metallicity regions reveals 
that about 76% of ordinary matter is in the form of hydrogen, about 24% (by mass) in '^He and trace abundances 
of ^He (- 10^^), deuterium (~ lO'^) and "^Li (- lO^^"), all relative to hvdrogenfl^ IT^. These elements were 
formed during the first three minutes of the Universe, while heavier elements (metals) are produced in the 
interior of stars and in astrophysical processes during supernovae explosions. 

• The concordance model: dark matter and dark energy. In the last few years there has been a wealth of 
observational evidence from CMB, large scale structure and high redshift supernovae Ia|3| data that leads to 
the remarkable conclusions that i) the spatial geometry of the Universe is fiat, ii) the Universe is accelerating 
today, and iii) most of the matter is in the form of dark matter. Current understanding of cosmology is based 
on the concordance or AC DM model in which the total energy densi ty o f the Universe has as main ingredients: 
5% of baryonic matter, 25% of dark matter and 70% of dark enerot/ [isl Hp. llTi ITsI] . The present observations 
indicate that the dark energy can be described by a cosmological constant [l9l|. 



B. The building blocks 

The main building blocks for a theory of the standard cosmology are: 

• Gravity: Classical general relativity provides a good description of the geometry of space time for scales 
I ^ Ipi ^ lO^'^^cm or time scales t ^ tpi ^ 10~''^s, or equivalently energy scales well below the Planck scale 
Mpi ~ lO^^Gev. A consistent quantum theory of gravity unified with matter describing the physics at the 
Planck scale is yet to emerge. 

Homogeneity and isotropy for a spatially flat space-time lead to the Friedmann-Robertson- Walker (FRW) metric 



ds^ = dt^ - a^it) dx^ (2) 

where t is the comoving time (the proper time of a comoving observer). Physical scales are stretched by the 
scale factor a{t) with respect to the comoving scales 

^physii) ~ ^(^) ^com ■ (3) 
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A physical wavelength redshifts proportional to the scale factor [eq.©], therefore its time derivative obeys the 
Hubble law lphys{t) = H{t) lphys{t) = lphys{t) I duit)- At equilibrium temperature decreases as the universe 
expands as 

In the homogeneous and isotropic FRW universe described by eq.©, the matter distribution must be homoge- 
neous and isotropic, with an energy momentum tensor with the fluid form 

{Ti:) = diag[p, -p, ^p, ~p] , (5) 

where p, p are the energy density and pressure, respectively. In such geometry the Einstein equations of general 
relativity reduce to the Friedmann equation, which determines the evolution of the scale factor from the energy 
density 



a{t) 



^H^t) = ^. (6) 



where Mpi = I/VSttG = 2.4 x 10^*^ GeV = 0.434 x 10"''' g. A spatially flat Universe has the critical density 

Pc = 3 Mil Hi = 1.88 /i^ lO^^^g/cm^ . (7) 

where Hq is the Hubble constant today. The energy momentum tensor conservation reduces to the single 
conservation equation, 

p + 3(p + p)-=0 (8) 
a 

The two equations 10 and © can be combined to yield the acceleration of the scale factor, 

a 1 



a 6 Mil 



(p + 3p) (9) 



which will prove useful later. In order to provide a close set of equations we must append an equation of state 
P = p{p) which is typically written in the form 

p = w{p) p (10) 
The following are important cosmological solutions: 



Cosmological Constant =^ ui = — 1 : de Sitter expansion p = constant ; a{t) cx e'^* ; H 

Radiation ^ w = — : Radiation domination ^ p oc a^^ ; a{t) oc (11) 
3 

Non — relativistic(cold) Matter w = : Matter domination => p oc ; a{t) oc 
Furthermore, we see from eqs.® and Ijlll) that accelerated expansion takes place if w < —1/3. 

• The Standard Model of Particle Physics: the current standard model of particle physics, experimentally 
tested with remarkable precision describes the theory of strong (QCD), weak and electromagnetic interactions 
(EW) as a gauge theory based on the group SU{3)c SU{2) (g) U{1)y- The particle content is: three generations 
of quarks and leptons: 

u\fc\ft\ ( Ve \ ( \ ( 



dj\sj\bj^\e J \ p J \ T 

vector Bosons: 8 gluons (massless) , Z°, VK± with masses Mz = 91.18 ±0.02 GeV and Mw = 80.4 ±0.06 GeV, 
the photon (massless) and the scalar Higgs, although the experimental evidence for the Higgs bosons is still 
inconclusive. 
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Current theoretical ideas supported by the renormahzation group running of the coupUngs propose that the strong, 
weak and electromagnetic interactions are unified in a grand unified theory (GUT) at the scale Mqut ^ 10^® Gev. 
Furthermore the ultimate scale at which Gravity is eventually unified with the rest of particle physics is the Planck 
scale Mpi ~ lO-'^^Gev. Although there are proposals for the total unification of forces within the context of string 
theories, their theoretical understanding as well as any experimental confirmation is still lacking. However, the physics 
of the standard model of the strong and electroweak interactions that describes phenomena at energy scales below 
~ 100 GeV is on solid experimental footing. 

The connection between the standard model of particle physics and early Universe cosmology is through Einstein's 
equations that couple the space-time geometry to the matter-energy content. As argued above, gravity can be studied 
semi-classically at energy scales well below the Planck scale. The standard model of particle physics is a quantum 
field theory, thus the space-time is classical but with sources that are quantum fields. Semiclassical gravity is defined 
by the Einstein's equations with the expectation value of the energy-momentum tensor T^'' as sources 

2^ Mil ^ ' 

The expectation value of T^^ is taken in a given quantum state (or density matrix) compatible with homogeneity and 
isotropy which must be translational and rotational invariant. Such state yields an expectation value of the energy 
momentum tensor with the fluid form eq.lOl. 

Through this identification the standard model of particle physics provides the sources for Einstein's equations. All 
of the elements are now in place to understand the evolution of the early Universe from the fundamental standard 
model. Einstein's equations determine the evolution of the scale factor, the standard model provides the energy 
momentum tensor and statistical mechanics provides the fundamental framework to describe the thermodynamics 
from the microscopic quantum field theory of the strong, electromagnetic and weak interactions. 

C. Energy scales, time scales and phase transitions 

Energy Scales: 

While a detailed description of early Universe cosmology is available in several books a broad- brush picture 

of the main cosmological epochs can be obtained by focusing on the energy scales of particle, nuclear and atomic 
physics. 

• Total Unification: Gravitational, strong and electroweak interactions are conjectured to become unified and 
described by a single quantum theory at the Planck scale ^ 10^^ Gev. There are currently many proposals that 
seek to provide such fundamental description such as string theories, however, their theoretical consistency is 
still being studied and experimental confirmation is not yet available. 

• Grand Unification: Strong and electroweak interactions (perhaps with supersymmetry) are expected to be- 
come unified at an energy scale ~ 10^^ Gev corresponding to a temperature T ~ 10^^ K under a larger gauge 
group G, for example S'?7(5), S'O(IO), -Eg, which breaks spontaneously G SU{3)c(^ SU{2)iS)U{1)y at a scale 
below grand unification. There are very compelling theoretical reasons for the existence of the GUT scale such 
as the merging of the running coupling constants of the strong, electromagnetic and weak interactions, shown 
in fig. ^for the minimal supersymmetric standard model (MSSM). Yet another reason is the explanation of 
the small neutrino masses via the see-saw mechanism in terms of the ratio between the weak and the grand 
unification scale. 

• Electroweak : Weak and electromagnetic interactions become unified in the electroweak theory based on the 
gauge group SU{2) (g) U{1)y- The weak interactions become short ranged after a symmetry breaking phase 
transition SU{2) ® U{1)y ~^ U{l)em at an energy scale of the order of the masses of the vector bosons, 
corresponding to a temperature Tew ^ 100 GeV ~ 10^^ K. At temperatures T > Tew the symmetry is restored 
and all vector bosons are (almost) massless (but for plasma effects that induce screening masses). For T < Tew 
the vector bosons that mediate the weak interactions (neutral and charged currents) , acquire masses 
through the Higgs mechanism while the photon remains massless, corresponding to the unbroken U{1) abelian 
symmetry of the electromagnetic interactions. Thus Tew determines the temperature scale of the electroweak 
phase transition in the early Universe and is the earliest phase transition that is predicted by the standard 
model of particle physics. The understanding of the symmetry breaking sector of the standard model is one of 
the primary goals of the Large Hadron Collider at CERN. The standard model has the necessary ingredients to 
explain the origin of the baryon asymmetry, leading to the possibility that the asymmetry between matter and 
antimatter was produced at the electroweak scale (see section lv)l . 



FIG. 1; Evolution of the weak ai, electromagnetic Q2 and strong as couplings with energy in the MSSM model. Notice that 
only 03 decreases with energy (asymptotic freedom) extracted from ref. |2ll| 

• QCD: The strong interactions have a typical energy scale Aqcd ~ 200 MeV at which the coupling constant 
becomes strong [as ~ C'(l)]- This energy scale corresponds to a temperature scale Tqcd ~ lO^^K. QCD 
is an asymptotically free theory, the coupling between quarks and gluons becomes smaller at large energies, 
but it diverges at the scale Aqcd. For energy scales below Aqcd QCD is a strongly interacting theory and 
quarks and gluons are bound into mesons and baryons. This phenomenon is interpreted in terms of a phase 
transition at an energy scale Aqcd or Tqcd- For T > Tqcd the relevant degrees of freedom are weakly 
interacting quarks and gluons, while below are hadrons. This is the quark-hadron or deconfinement-confinement 
phase transition. In the limit of massless up and down quarks, QCD features an SU{2)]^ (g) SU{2)^ chiral 
symmetry, which is spontaneously broken at about the same temperature scale as the confinement-deconfiment 
transition. Pions are the (quasi) Goldstone bosons emerging from the breakdown of the chiral symmetry 
SU{2)]^ (g) SU{2)fi SU [2)11^1^. The QCD phase transition(s) are the last phase transition predicted by 
the standard model of particle physics. The high temperature phase above Tqcd, with almost free quarks and 
gluons (because the coupling is small by asymptotic freedom) is a quark-gluon plasma or QGP. Experimental 
programs at CERN (SPS-LHC) and BNL (AGS-RHIC) study the QCD phase transition via ultrarelativistic 
heavy ion collisions (URHIC) and a systematic analysis of the data gathered at SPS and RHIC during the last 
decade has given an optimistic perspective of the existence of the QGP 22, 23] (see section 4 below). 

• Nuclear Physics: Low energy scales that are relevant for cosmology are determined by the binding energy 
of light elements, in particular deuterium, whose binding energy is ^ 2 Mev corresponding to a temperature 
T ~ 10^° K. This is the energy scale that determines the onset of primordial nucleosynthesis. The first step in the 
network of nuclear reactions that yield the primordial elements is the formation of the deuteron via n+p <-> (i+7. 
The large number of photons per baryon results in that high energy photons in the blackbody tail dissociate the 
deuterons formed in the forward reaction until the temperature becomes of the order of Tms ~ 0.1 MeV. Once 
deuterons are formed, a network of nuclear reactions results in that all neutrons end up in nuclei, mainly helium, 
resulting in a helium abundance of about 25%|0,0|. The nature of the nuclear forces suggests the possibility 
of a liquid-gas phase transition in cold nuclear matter at an energy ~ 15 — 20 MeV discussed in section IVIII El 

• Atomic Physics: A further very important low energy scale relevant for cosmology corresponds to the binding 
energy of hydrogen ~ 10 eV. This is the energy scale at which free protons and electrons combine into neutral 
hydrogen or 'recombination'. The large number of photons per baryon results in that recombination actually 
takes place at an energy scale of order 0.3 eV, at about 400000 years after the beginning of the Universe. At 
this time when neutral hydrogen is formed the Universe becomes transparent. This event determines the last 
scattering surface, after neutral hydrogen is formed photons no longer scatter and travel freely. These are the 
photons measured by CMB experiments today. 

Time Scales : An important ingredient of modern standard cosmology is a brief but explosive early period of 
inflation during which the scale factor grows exponentially as a{t) = e [see eq. (g it]- W MAP[il[ll[ia[lg yields 
an upper bound on the energy scale of the inflation [for a detailed discussion see sec. IIV C| H < 10^^ GeV. In order to 
solve the entropy and horizon problems, the inflationary stage must last a time interval St so that H St ^ 60, hence 
the inflationary stage lasts a time scale 



Aim/ ~ 10 sees 



(13) 
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Field models of inflation are discussed in sec. IIV CI The inflationary stage is followed by a radiation dominated era 
(standard hot big bang) after a short period of reheating dm'ing which the energy stored in the field that drives inflation 
decays into quanta of many other fields, which through scattering processes reach a state of local thermodynamic 
equilibrium. 

Once local thermodynamic equilibrium is reached, a very detailed picture of the thermal history of the Universe 
emerges combining statistical mechanics with the basic ingredients described above during the first ~ 10000 

years of the Universe and after the infiationary stage that lasted ~ 10~'^'*secs, the Universe was radiation dominated 
expanding and cooling (almost) adiabatically. As a consequence the entropy S cx V{t) T^{t) oc Vq [a{t) T{t)]'^ is 
almost constant according to eq.Q) but for the change in the number of relativistic degrees of freedom. Therefore, 
for a relativistic equation of state p = p/3 with p being the pressure and p the energy density, p{t) oc T'^{t). Such 
equation of state yields the following evolution of the scale factor as a function of time 

where g is the number of relativistic degrees of freedom, which is also a function of temperature 10 < 5 < 100 for 
0.5 MeV ^ T < 300 GeV. The above expression yields a simple dictionary that allows to translate temperature (or 
energy scale) into time scales, namely 

MeV 10i°K 

(15) 



y/t{sec) ^t(sec) 

This simple dictionary allows to establish the time scales at which the standard model of particle physics predicts 
phase transitions as well as an estimate of whether the transitions are likely to occur in LTE or not. The electroweak 
transition would have ocurred at T '--^ 100 GeV at a time scale tEw ~ 10~^^secs and the QCD phase transition at 
T - 170 MeV at tqcD ^ 10"^ sees. 

Local Thermal Equilibrium (LTE) or Nonequilibrium: 

Whether a phase transition occurs in or out of local thermodynamic equilibrium (LTE) depends on the comparison 
of two time scales: the cooling rate and the rate of equilibration. In the early Universe or in ultrarelativistic heavy 
ion collisions, the rate of change of temperature is determined by the expansion rate of the fiuid. The rate of cooling 
by cosmological expansion follows from eq.Q T{t)/T{t) = ~a{t)/a{t) — —Hit). Collisions as well as non-collisional 
processes contribute to establish equilibrium with a rate T. Local thermodynamic equilibrium ensues when T > H{t), 
in which case the evolution is adiabatic in the sense that the thermodynamic functions depend slowly on time through 
the temperature. When the cosmological expansion is too fast, namely H(t) 3> T, local thermodynamic equilibrium 
cannot ensue, the temperature drops too fast for the system to have time to relax to LTE and the phase transition 
occurs via a quench from the high into the low temperature phase. 

While a detailed understanding of the relaxational dynamics requires an analysis via quantum Boltzmann equations, 
a simple order of magnitude estimate for a coUisional rate is given by the ensemble average F ~ {a nv), where a is 
a scattering cross section, n is the density of scatterers and v the average velocity. For electromagnetic scattering a 
typical cross section is of order ~ cx^/Q^ with the transferred momentum, at high temperature single photon 
exchange yields the estimate (Tem ~ ct'^/T'^, the density of ultrarelativistic degrees of freedom n ^ and v ^ 1 
yields Tem ^ cP' T. In QCD, simple gluon exchange yields the estimate Tqcd ^ al T. Comparing with H given 
by ea. H14|) . it is found that the strong interactions are in LTE for T < 10^^ GeV and electromagnetic interactions 
are in LTE for T < 10^^ GeV. A similar estimate is obtained for the weak interactions: a typical scattering process 
with energy transferred E ^ AI\y has a scattering cross section a ^ G\ whereas if E » Mw, o ^ 9^ 1^"^ i 
therefore in a thermal medium with E ^ T and density of relativistic particles n a. typical weak interaction 

reaction rate is Vew ^ T iax T ]^Iw,z, and Tew G% for T <^ Mw,z- In this latter temperature regime 
the ratio Tew/H ^ (T/MeV)'^, hence the weak interactions fall out of LTE for T < 1 MeV. This simple analysis, 
while providing an intuitive order of magnitude estimate for the relaxation time scales, neglects several subtle but 
important aspects that must be studied on a case-by-case basis: 

• Screening and infrared phenomena: the estimate for the relaxational rates F invoked the exchange of a vector 
boson or relativistic degrees of freedom. In a medium at high temperature and or density there are important 
screening effects and infrared phenomena that change these assessments both quantitatively and qualitatively 
and depend on whether the gauge symmetry is abelian or not|24j. 

• Critical slowing down: condensed matter experiments reveal that systems that undergo second order (or in gen- 
eral continuous) transitions exhibit a slowing down of relaxational dynamics of long wavelength fiuctuations near 
the critical point of a phase transition. Such is the case in ferromagnets, superconductors and superfluids. While 
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short wavelength fluctuations remain in LTE through the transition through microscopic scattering mechanisms, 
long wavelength fluctuations feature slower relaxational dynamics and even cooling rates far smaller than micro- 
scopic relaxational rates produce quenched transitions and departures from equilibrium. Critical slowing down 
in classical models of critical phenomena is fairly well understood '2 5', '26', '2^ , but a similar level of understanding 
in quantum field theories at extreme temperature and density is now emerging _28J. 

• Strong first order transitions and metastable states: In a strong first order transition (see below) the system is 
trapped in a metastable state which is a local minimum of the free energy but not the global one. Within this 
local, metastable minimum, collisions can bring the system to LTE, but the metastable state will eventually 
decay and fall out of LTE by the non-perturbative process of nucleation, which is described below in detail. 



D. Phase transitions: early Universe vs. accelerator experiments 

The control variables in a collision experiment are the beam energy and the luminosity. For a phase transition to 
be achieved in an accelerator experiment an environment with a temperature close to the transition value Tc must 
be formed in the collision region. The blackbody relation, valid for ultrarelativistic particles in equilibrium yields the 
energy density-temperatwce relation 

e = CT'^, (16) 

where the constant C depends on the number of degrees of freedom. For the electroweak phase transition with 
Tc ~ 100 GeV, the energy density that must be deposited in the collision region to achieve the conditions for a phase 
transition is e '--^ 10^° GeV/fm'^. This energy density is about 10^^ times larger than that of nuclear matter. For the 
QCD phase transition with Tc ~ 177 MeV, e '--^ GeV/fm'^ which is achieved in ultrarelativistic heavy ion collisions 
at SPS-CERN and RHIC-BNL. Therefore, from the perspective of studying particle physics phase transitions with 
accelerator experiments the only realistic possibility for the foreseeable future is the QCD transition with URHIC at 
RHIC and the forthcoming LHC. Hence the potential observables from phase transitions in the early Universe before 
the QCD scale must be inferred indirectly from the aftermath. In section rVIII Dl we compare the conditions for the 
QCD phase transition both in the early Universe and in URHIC to assess whether current experiments reproduce the 
conditions that prevailed about 10 /i sees after the beginning of the Universe. 



III. PHASE TRANSITIONS: EQUILIBRIUM AND NON-EQUILIBRIUM ASPECTS: 
A. Equilibrium aspects: free energy, effective potentials and critical phenomena. 

Phase transitions are broadly characterized as either second or first order pM l2^ l27j . In the Landau theory of phase 
transitions the order parameter plays a central role. In particle physics the order parameter is the expectation value 
of a spin zero field in the state that extremizes the free energy. A non-zero expectation value for a non-zero spin field 
entails the breakdown of rotational symmetry. A particularly illuminating example is that of ferromagnetic materials 
where the order parameter is the total magnetization of the sample. For these materials the magnetization vanishes 
above the Curie temperature while a nonzero spontaneous magnetization emerges below such critical temperature. 
In the standard model, the order parameter is the expectation value of the neutral component of the Higgs {(fio), in 
QCD the chiral order parameter is the expectation value of the pseudoscalar density (tp^^ip). For the confinement- 
deconfincment phase transition lattice studies show a rapid variation of the free energy as a function of temperature 
and the Polyakov loop may play the role of order parameter (see below). 

At zero temperature and chemical potential the state of lowest free energy is the vacuum state of the theory. An 
important concept to study the nature of the phase transition is that of the effective potential or Gibbs free energy 
which is a function of the expectation value of the scalar field. 

Consider adding to the Hamiltonian of the theory H a term J ip where J is a constant and (f is the scalar field 
whose expectation value is the order parameter. The total Hamiltonian is H[J] — H + J (p. The partition function is 
given by 

Z[J,T] = Tre-^ (17) 

introducing the Helmholtz free energy 



W[T, J] = -T \nZ[J,T] , 



(18) 
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the order parameter $ is given by 



5W[T, J] 
5J 



(19) 



inverting this relation one finds J(<i>). The Gibbs free energy or effective potential V[^] follows from a Legendre 
transform 



= VF[J($)] - J($)$ 



(20) 



and it is a function of the order parameter and other intensive thermodynamic variables such as temperature, chemical 
potential, etc. It is a very powerful concept that provides information on the equilibrium thermodynamic aspects and 
the phase structure of the theory and the possible transitions between them. 




FIG. 2: Effective potentials for a second order phase transition with the breakdown of a discrete symmetry (left panel) and 
for a first order phase transition (right panel). The filled circles in the right panel figure display the equilibrium minima for 
T > Tc and T — Tc and the metastable minimum for T < Tc- 



Fig. 121 displays the typical effective potentials for a second order (left panel) and first order (right panel) phase 
transitions. In the case of a second order transition, the order parameter in the state of minimum free energy vanishes 
continuously as the temperature approaches the critical. In second order phase transitions the second derivative of the 
free energy with respect to temperature is divergent at the critical temperature. Near the critical temperature the free 
energy per unit volume, or equivalently the pressure varies as p{T) ~ \T — Tc\'^^" with a the thermal critical exponent. 
The order parameter itself vanishes as \T — Tc\^ for T ^ Tc from below, with (3 > and response functions feature 
power-law singularities. In contrast to this case, in a first order phase transition the value of the order parameter 
and the position of the minimum of the free energy jump discontinuously at the critical temperature. In fig. |21the 
black circles in the right panel show the behavior of the order parameter as the temperature is lowered from above 
to below Tc- The state that minimizes the effective potential changes suddenly from the right to the left well as the 
temperature is lowered below Tc. For T > Tc the right well in the right panel remains a local minimum of higher 
free energy, thus a metastable state. The difference in free energies between the local and the global minima yields a 
latent heat which is released upon the decay of the metastable state. 

The effective potential also provides information on the excitation spectrum. Consider a state characterized by a 
given value of the order parameter $, the frequency of harmonic excitations of wavevector k around this state is given 
by 

nik;^)^Vk^ + V"i^) , V"i^) = ^^^. (21) 

Most of the phase transitions in particle physics models involve the spontaneous breakdown of a global symmetry 
and the order parameter transforms covariantly under this symmetry. An expectation value of the order parameter 
in the state of lowest free energy indicates the spontaneous breakdown of the symmetry. Different phase transitions 
in particle physics involve different order parameters: for the electroweak transition the order parameter is the 
expectation value of the Higgs doublet along the neutral component, = {ifo}- This expectation value breaks the 
symmetry SU{2) (g) U{1)y U{l)em, the three Goldstone bosons emerging from the global broken symmetry make 
up the longitudinal components of the W^, vector bosons via the Higgs mechanism. In QCD with massless quarks 
the expectation value of the chiral density (^'7^^') with ^E* an up-down quark doublet breaks the SU{2)l (8) SU{2)fi 
symmetry of the QCD Lagrangian down to 5L/(2)l+_r, with pions emerging as the triplet of Goldstone bosons. 
Because up and down quarks are not exactly massless, pions are pseudo Goldstone bosons. The QCD confinement- 
deconfinement transition does not have an obvious order parameter. However, in absence of dynamical quarks, for 
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example in an SU{N) Yang- Mills theory of gluons, the Polyakov loop 

L{x) = ^ ( TrT'exp (ig Ao(f,T) dr ) \ (22) 



N 







has been shown to be a suitable gauge invariant order parameter. In the above expression g is the gauge coupling, T is 
the temperature Aq is the temporal component of the SU{N) gauge potential, and V is the path ordering symbol, the 
line integral is carried out in Euclidean time. The Polyakov loop is usually interpreted as the free energy of an infinitely 
heavy test quark. It vanishes in the confined phase for T < Tc and is non- vanishing in the deconfined phase T > Tc- 
An effective potential in terms of the Polyakov line has recently been proposed to study the confining-deconfining 
phase transition j22i|. When light dynamical quarks are included there is no longer an obvious order parameter for the 
confinement-deconfinement phase transition because light dynamical quarks screen the potential from gluon exchange. 

Second order phase transitions in equilibrium and critical phenomena: Consider a continuous (second 
order) phase transition _25j described by an effective potential akin to the left panel in fig. [5] in the case in which 
the temperature falls from a value larger than the critical value Tc. If the cooling rate is much smaller than the 
relaxation rate the evolution will be in LTE, but as the temperature approaches the critical, collective long-wavelength 
fluctuations develop and become strongly correlated. Near Tc these collective long-wavelength fluctuations become 
massless (critical), their correlation function becomes scale invariant and large regions behave coherently with strong 
correlations over arbitrarily large scales. This is the hallmark of critical phenomena at second order phase transitions 
which are characterized by response functions that feature singularities in \T — Tc\ in terms of power laws with 
critical exponents. The correlation function of the order parameter field $ is given by ($(f)$(0)) cx e-^/«(^). The 
correlation length ^(T) diverges at the critical temperature as ^(T) ^ \T — Tc\~'^ ^ with v > Q being a critical exponent. 
At this point the system becomes correlated over large distances and the spectrum of long-wavelength fluctuations 
becomes scale invariant. This is a very important aspect of critical phenomena associated with second order phase 
transitions. Near the critical point where the correlation length diverges the long distance physics is universal and the 
static aspects can be described in terms of a Landau-Ginzburg low energy effective theory. The different universality 
classes that characterize the different critical exponents, are determined by few properties of the system such as the 
dimensionality, the symmetry of the order parameter, the number of independent fields, etc. The description of static 
critical phenomena in terms of the Landau-Ginzburg approach has been confirmed by a wealth of experiments in 
condensed matter systems. A relevant example of a Landau-Ginzburg description for the case of a scalar field, whose 
expectation value breaks a discrete symmetry ip ^ —ip is described by the finite temperature effective potential given 
by 

V[<P; T] = I (T - Te) $2 + ^ $4 (-23) 

and featured in fig|21 The concept of universality classes that describes critical phenomena for a wide variety of 
systems is of fundamental importance. Relevant to the discussion below is the case of the chiral phase transition 
in QCD, which for massless up and down quarks, is described by the same universality class as the Heisenberg 
ferromagnet with 0(4) symmetry [30l|. 

An important dynamical aspect of second order phase transitions is that near the critical region the relaxation time 
scale of long- wavelength fluctuations feature critical slowing down psl l26l l27| . Critical slowing down has been studied 
theoretically and experimentally in condensed matter systems, where a large body of experimental work confirms the 
critical slowing down of these fluctuations near second order phase transitions in classical systems in agreement with 
theory. The study of critical slowing down in quantum field theory has recently began to receive attention ^2^ . For 
any finite cooling rate, long-wavelength modes will be quenched through the phase transition. The effective potential 
provides only an equilibrium description but does not address the issue of the dynamics of the transitions between 
different phases. 



B. Non-equilibrium aspects: spinodal decomposition and nucleation. 

As discussed in sec lll Cl whether a phase transition occurs in LTE or not depends on whether, T > H{t) or F < H{t), 
respectively. In the first case the transition occurs adiabatically (similar to the Minkowski case) while in the latter 
case, the phase transition occurs via a quench from the high into the low temperature phase. The dynamics are 
different depending on whether the equilibrium effective potential describes a second or first order transition. 
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1. Second order case: spinodal decomposition 



Let us first consider tlie simpler case of a symraetry brealcing second order transition wlicrc the symmetry tliat breaks 
spontaneously is discrete. To focus on a simple, yet relevant example, consider the form for the finite temperature 
effective potential given by ea. (|23|l . For T > the state of minimum free energy corresponds to $ = 0, whereas for 
T < Tc the state of minimum free energy corresponds to either the right or left equilibrium minima in fig. [2 namely 



<f> = ±<i>p with $p 



± 



^(1 



^) . The inflection points at which V"{^,T) are at <i> = ±^s with 



^^(1 — ^) . The states with — <f>s < ^ < +$s are thermodynamically unstable. In order to understand the non- 
equilibrium aspects of a rapid phase transition it is convenient to plot the coexistence lines at which thermodynamic 
equilibrium states coexist and the spinodal line that limits the region of thermodynamic stability in the T — (p 
plane. These are shown in fig. |21 The coexistence line is given by the relation between $ and T for the minima 
of the potential, namely Tcoex{'^) = ^^/a + Tc, and similarly the spinodal line is determined by the condition 
F"($,T) = 0, namely the inflexion points of the effective potential, which yields Tsp{<^) — —b $^/(3 a) +Tc. Fig. 
13 depicts a quench from a high temperature state with vanishing order parameter to a low temperature state when 
the equilibrium state corresponds to a broken symmetry. In classical statistical mechanics out of equilibrium the 
dynamics of the evolution out of equilibrium is studied via purely dissipative equations of motion, the Cahn-Hilliard 
equations |25i. i26i. i27j . The initial stages of spinodal decomposition in classical systems are characterized by growth 
of unstable long- wavelength perturbations which lead to the formation of domains and coarsening j2^ |23|- Beautiful 
experiments in phase separation and spinodal decomposition in binary fluids have provided a spectacular confirmation 
of this mechanism as well as the dynamical aspects predicted by the Cahn-Hilliard theorvlsij. In quantum field 
theory the dynamics is determined by the unitary time evolution of an initially quantum state or density matrix. 
Quantum spinodal decomposition is studied in ref. [s^ . Is^ and the non- linear backreaction of unstable fluctuations 
was studied in by implementing a non-perturbative self-consistent method. States between the spinodal lines 
are thermodynamically unstable to small amplitude long- wavelength perturbations, since < in the spinodal 

region +*i's, the frequency of small amplitude fluctuations around this state become imaginary for 

wavevectors k < \V"{^)\, while fluctuations with k > \V"{^)\ remain stable. Fluctuations with wavevectors in the 
spinodal band k < y/\V"{0)\ grow as e'^'' * with jk = — k'^- These unstable modes lead to the formation of 

domains whose size is determined by a time dependent correlation length ^{t). Inside these domains the expectation 
value attains the equilibrium values ±$e at the final temperature. The formation of correlated domains of each 
phase is the hallmark of the process of phase separation via spinodal decomposition. The dynamics of this process is 
non-perturbative: initially long wavelength fluctuations grow exponentially until they become non-linear and react 
back on the dynamics shutting off the instabilities. In scalar field theories it is found that during the initial stages 
of phase separation f (i) cx t. The above discussion focused on the case of discrete symmetry that is spontaneously 
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FIG. 3: A quench into the spinodal region. 



broken. When the broken symmetry is continuous a wealth of fascinating non-perturbative topological excitations 
emerge, such as vortices, monopoles, and textures j2J]- 
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2. First order phase transitions: nucleation 



The dynamics of first order phase transition is different: in this case the thermodynamic state of the system is 
a trapped metastable siaiejlocally stable under small thermodynamic perturbations, but with a higher free energy 
than the global minimum |2a. l27l| . This is the situation of the liquid-gas phase transition for example in water. The 
metastable state is stable under small amplitude perturbations (oscillations in the right-most well in fig. ^ but decays 
via spontaneous large amplitude fluctuations corresponding to bubbles of the stable phase immersed in a host of the 
metastable state. Consider a spherical such bubble or radius R, where inside of the bubble the value of the order 
parameter is that of the globally stable state, but outside is that of the metastable state. Since the globally stable 
state has a lower free energy than the metastable state, there is a gain in the free energy given by —AttR^AV/^ 
where is the difference in the effective potential between the stable and the metastable state. Because the order 
parameter is inhomogeneous in this configuration there is an elastic contribution to the free energy from the gradients 
of the order parameter field which is proportional to the surface of the bubble, because this is the region in which the 
spatial derivatives of the order parameter are non- vanishing. This elastic contribution is positive and given by AnE? a 
where a is the surface tension, thus the total change in the free energy for such inhomogeneous bubble configuration 
is 

E[R] = Attract - — R^AV . (24) 
o 

This change in the free energy is depicted in fig. 0J 
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FIG. 4: E{R) vs. R, arbitrary units. 



Bubbles with radii R < R* — 2a/ AV shrink while those with R > R* — 2(t/AV grow and convert the metastable 
phase into the globally stable phase releasing the latent heat. The phase transition completes via percolation of these 
growing bubbles. The probability of a bubble to appear spontaneously in the heat bath in equilibrium at temperature 
T is given by 

167r fj"^ 

P{R*)=De'^ ■ E[R*] = —j^ (25) 

where the prefactor D depends on the spectrum of small amplitude fluctuations around the critical bubble configura- 
tion. The rate of decay per unit volume and per unit time of the metastable state is given bv p5ll2^l27| 

E[R*] 

r = Toe-— (26) 



where Fq depends on the prefactor D as well as in general on transport properties such as the heat conductivity [2y,|23|. 
These processes which require large amplitude fluctuations to overcome a potential barrier, which in this case is 
determined by are called thermally activated and the decay rate features the Boltzmann suppression factor 

corresponding to the energy of the configuration at the top of the barrier. Since the metastable state has a higher 
free energy than the globally stable state, the nucleation process that results in the conversion of the metastable into 
the stable phase must release the latent heat. 



3. Liquid-gas phase transition: phase coexistence 



Another important and very relevant example of a first order phase transition is the liquid-gas transition. This 
transition is exemplified by an equation of state of the Van der Walls form and results in general when the microscopic 
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interactions have a short range repulsive and a long range attractive components. Such is the case of nuclear forces 
between nucleons (see section IVlll Jil|) . Fig. |S1 displays a few isotherms for a typical equation of state in terms of 
pressure P vs. density p for a liquid-gas type transition. For T > Tc the isotherms are single valued, there is 
only one value of density p for a given pressure P and the system is in one phase. However, for T < Tc, for a 
given value of the pressure there are three values of the density. The derivative dP/dp\T = c^t^ where Cs.t is the 
isothermal speed of sound for hydrodynamic modes. Obviously, we see that < in the region between pa and 
Pb in figEl This corresponds to an imaginary speed of sound, therefore to unstable density fluctuations that grow 
exponentially with time as e*^ I'^^.tI Therefore, such unstable thermodynamic states and the value of the pressure in 
this region do not correspond to a thcrmodynamically stable equilibrium state. The two different values pa,b describe 
two different phases: the phase corresponding to pA features a small speed of sound, hence a highly compressible 
medium, this is the gas phase, while the phase described by pb features a large speed of sound, hence a rather small 
compressibility, namely a liquid. The dashed line at constant pressure in fig. |S1 describes the coexistence of these 
phases in thermodynamic equilibrium at the same pressure, temperature and chemical potentials. The position of this 
line is determined by the equal area law or Maxwell construction. Along the coexistence curve the state of the system 
is a heterogeneous mixture of both phases with a composition given by the lever rule: a state with global density p 
is a mixed phase with islands of gas with density pA in a host of liquid with density ps with p = C,pA + (1 — Orb 
where C — {pA — Pb) / {p — Pb) < 1 is the proportion of the phases. A change of the global density in this region 
results in a conversion of part of one of the phases into the other at constant pressure and temperature. Therefore 
in the coexistence region or mixed phase small variations in the density do not result in a change in the pressure or 
temperature, just in the isothermal and isobaric conversion of phases, hence the isothermal speed of sound vanishes 
along the coexistence curve or mixed phase. This flat portion of the isotherm in the mixed phase is a soft part of 
the equation of state since the isothermal speed of sound vanishes. This property of the mixed phase or coexistence 
region for first order phase transitions has important consequences both in cosmology as well as in ultrarelativistic 
heavy ion collisions to be discussed below. 

The region pa < p < Pb corresponding to dP/dp < is the spinodal region: small pressure (or density) perturbations 
of a homogeneous phase are unstable and grow. The regions pA < p < pa ', Pb < p < Pb are actually metastable. 
Consider the following experiments: a) begin at T > with a state in thermodynamic equilibrium with a single 
phase at density Pa < P < Pb^ and quench the system in temperature by lowering the temperature to T < Tc on 
a very short time, b) repeat the experiment but now quenching the system from T > Tc to T < Tc at constant 
density pb < p < pb, what is the dynamics in both cases? In the case (a) the initial state is a single (gas) phase, 
homogeneous of density p which is quenched inside the spinodal region where small amplitude density perturbations 
of a homogeneous phase arc unstable. The spinodal instabilities lead to a fragmentation of the homogeneous phase 
and the formation of domains of the liquid phase, these domains grow and separate, leading to the process of phase 
separation and an heterogeneous mixed phase in coexistence with proportions of phases given by the lever rule. In 
the case (b) the initial homogeneous state (gas) with global density pb < p < pb quenched into a region in which 
the homogeneous phase is mechanically stable since dP/dp > 0, however at this temperature the free energy of the 
liquid is smaller, hence the state is metastable. It will decay to a mixed state of gas and liquid along the coexistence 
line by nucleation: a large amplitude fiuctuation corresponding to a droplet larger than the critical, the growth and 
percolation of these droplets eventually leads to the mixed phase on the coexistence curve. Spontaneous nucleation 
of a liquid droplet in the homogenous gaseous host is referred to as homogeneous nucleation, but also the presence of 
impurities can act as nucleation seeds, in which case the situation is referred as heterogeneous nucleation. 

This is in fact the principle of the cloud chamber, one of the first particle detectors where a gas was supercooled 
and cosmic rays (or other particles) passing through the gas induced the nucleation of droplets, namely these particles 
triggered the heterogeneous nucleation of the metastable phase. 




Since this state is continuously connected to the low density region, this single phase describes a gas. 
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FIG. 5: Left panel: equation of state for a liquid gas transition. The dashed line is the coexistence line (Maxwell construction), 
Pa, 6 determine the position of the spinodal line. Right panel: coexistence and spinodal lines in the T — p plane. 



The spinodal and coexistence regions end at the critical point, and for T > the system can only be in a single 
thermodynamically stable phase. Any path taking the system from a liquid into the gas phase above Tc is continuous. 
The critical point itself is of second order. If the system is brought near the critical temperature from above but 
always in LTE, the system develops large density fluctuations with a correlation length that diverges. These strong 
fluctuations in liquid-gas systems gives rise to the phenomenon of critical opalescence which is observed in the familiar 
example of water near its boiling point when bubbles of vapor are beginning to be formed and the liquid becomes 
opaque as a result of light scattering from the large vapor domains. 

A liquid-gas type transition is ubiquituous in systems that feature microscopic short range repulsion and a long range 
attraction. Such is indeed the case of nucleons in nuclei, semi-phenomenological effective nucleon interactions yield 
equations of state of the Van der Walls type and predict that nuclear matter undergoes a liquid-gas phase transition at 
a critical temperature ~ 15 — 20 MeV. The theoretical and experimental aspects of such phase transition in nuclear 
matter are discussed further in section IVIII El // QCD has a first order confinement-deconfinement (quark hadron) 
phase transition the soft-part of the equation of state corresponding to a mixed phase gives rise to an anomalously 
small speed of sound with important consequences both in the early Universe as well as in URHIC. A mixed phase 
during the quark-hadron transition in the early Universe would result in an almost vanishing speed of sound, as a 
result there are no pressure gradients that can hydrostatically balance the gravitational pull and gravitational collapse 
would ensue, perhaps resulting in the formation of primordial black holes. In URHIC, if the quark-hadron gas reaches 
coexistence as a mixed phase in equilibrium, again the small speed of sound (soft point) hinders the pressure gradients 
that drive the expansion and hydrodynamic flow, with potential experimental observables. 

Crossover: An alternative to a second or first order phase transition is a simple crossover in thermodynamic 
behavior without discontinuities or singularities in the free energy or any of its derivatives. If the crossover is 
smooth, then no out-of-equilibrium aspects are expected as the system will evolve in LTE. However, if the crossover 
is relatively sharp the situation may not be too different from a phase transition. There is now evidence suggesting 
that the standard model does not feature a sharp electroweak phase transition (either first or second order) but is a 
smooth crossover. Furthermore, current lattice studies indicate that for two light quarks (u,d) and one heavier quark 
(s) with a mass of the order of the QCD scale the confinement-deconfinement transition maybe a sharp crossover 
rather than cither a first or second order transition. These possibilities are discussed further below. 

IV. INFLATION AND WMAP 
A. Inflationary dynamics 

Inflation was originally proposed to solve several outstanding problems of the standard big bang model [s^ thus 
becoming an important paradigm in cosmology. At the same time, it provides a natural mechanism for the gen- 
eration of scalar density fluctuations that seed large scale structure, thus explaining the origin of the temperature 
anisotropies in the cosmic microwave background (CMB), as well as that of tensor perturbations (primordial gravi- 
tational waves). Inflation is the statement that the cosmological scale factor a{t) in eq.ll^J has a positive acceleration, 
namely, a{t)/a{t) > 0. Hence, eq.Q requires the equation of state w — p/p < —1/3. 

Inflation gives rise to a remarkable phenomenon: physical wavelengths grow faster than the size of the Hubble 
radius dn — a{t)/a{t) = 1/H{t), indeed 

^^^=Hit)^^ + d,''^, (27) 

Ea. (|27|l states that during inflation physical wavelengths become larger than the Hubble radius. Once a physical 
wavelength becomes larger than the Hubble radius, it is causally disconnected from physical processes. The inflationary 
era is followed by the radiation dominated and matter dominated stages where the acceleration of the scale factor 
becomes negative since p/ p = 1/3 in a radiation dominated era and p = in a matter dominated era [see eq. . With a 
negative acceleration of the scale factor, the Hubble radius grows faster than the scale factor, and wavelengths that were 
outside, can now re-enter the Hubble radius. This is depicted in figEl This is the main concept behind the infiationary 
paradigm for the generation of temperature fluctuations as well as for providing the seeds for large scale structure 
formation: quantum fluctuations generated early in the inflationary stage exit the Hubble radius during inflation, and 
eventually re-enter during the matter dominated era. The basic mechanism for generation of temperature anisotropies 
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FIG. 6: Logarithm of physical scales vs. logarithm of the scale factor. The causal horizon dn is shown for the inflationary 
(De Sitter), radiation dominated and matter dominated stages. The physical wavelengths for today's Hubble radius dn (today), 
and a galactic scale Xgai axe shown. 



as well as primordial gravitational waves through inflation is the following [l|-[7|: the energy momentum tensor is split 
into a the fluid component T^i^^^ [eq.jSjl] that drives the classical FRW metric plus small quantum fluctuations, 
namely T^^ = "^/hiid ^ quantum fluctuation of the matter fields induce a quantum fluctuation in the metric 

(geometry) 6 G'^'^ — 6T'^'^ /Mpi. In the linearized approximation the different wavelengths of the perturbations evolve 
independently. After a given wavelength exits the Hubble radius, the corresponding perturbation becomes causally 
disconnected from microphysical processes. Perturbations that re-enter the Hubble radius at the time of photon 
decoupling, about 400000 years after the beginning of the Universe, induce small fluctuations in the space time metric 
which induce fluctuations in the matter distribution driving acoustic oscillations in the photon-baryon fluid. At the 
last scattering surface, when photons decouple from the plasma these oscillations are imprinted in the power spectrum 
of the tem pera ture anisotropies of the CMB and seed the inhomogeneities which generate structure upon gravitational 
collapse p. l36l|. The horizon problem, namely why the temperature of the CMB is nearly homogeneous and isotropic 
(to one part in 10^) is solved by an inflationary epoch because the wavelengths corresponding to the Hubble radius 
at the time of recombination were inside the Hubble radius hence in causal contact during inflation. This mechanism 
is depicted in fig. While there is a great diversity of inflationary models, they generically predict a gaussian and 
nearly scale invariant spectrum of (mostly) adiabatic scalar (curvature) and tensor (gravitational waves) primordial 
fluctuations. These generic predictions of inflationary models make the infiationary paradigm robust. The gaussian, 
adiabatic and nearly scale invariant spectrum of primordial fiuctuations provide an excellent fit to the highly precise 
wealth of data provided by the Wilkinson Microwave Anisotropy Probe ( WMAP) [H [H [ll[l3 . WMAP has also 
provided perhaps the most striking validation of inflation as a mechanism for generating superhorizon fluctuations, 
through the measurement of an anticorrelation peak in the temperature-polarization (TE) angular power spectrum 
at ^ ~ 150 corresponding to superhorizon scales p^ll7l[l8| . 



B. Inflation and scalar field dynamics 



A simple implementation of the infiationary scenario is based on a single scalar field, the inflaton with a Lagrangian 
density 



C = a^(t) 



(v^)2 



2a^(t) 



V(ip) 



(28) 



pressui 

infiaton (p(t) field in the early universe [see sec. Ill Ij| are given by 



with V(ip) the inflationary potential. The energy density and the pressure for a spatially homogeneous and isotropic 

(29) 



• 2 '2 

p^\ + v(ip) , p = ^-y(^) 
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Inflation should last at least > 60 efolds in order to solve the entropy and horizon problems. This entails a 
slow evolution and small temporal derivatives for the inflaton (slow roll), namely ip^ <^ '^{f)- This implies p ~ 
—p ~ V{(p) ~ constant as the equation of state leading to a de Sitter universe with scale factor a{t) — e^*, H = 
^JV{ip)/[3 M'pi] [see eq. Jll ll]. This situation is achieved via a variety of inflationary scenarios, old, new, chaotic, 
hybrid inflation etc. (see [3|-[3 for discussions on different models). 

While inflationary dynamics is typically studied in terms of a classical homogeneous inflaton fleld, such classical fleld 
must be understood as the expectation value of a quantum field in an isotropic and homogeneous quantum state. In 
ref.[37, 38] the quantum dynamics of inflation was studied for inflation potentials which features a discrete symmetry 
breaking, V{lp) = — m^(y5^/2 + Xip'^ (new inflation) as well as unbroken symmetry potentials V{lp) — -Vrn^ip^ /2 + Xip'^. 
The initial quantum state was taken to be a gaussian wave function(al) with vanishing or non- vanishing expectation 
value of the field. This state evolves in time with the full infiationary potential which features a spinodal region for 
(p^ < w? /12\ in the broken symmetric case. Just as in the case of Minkowski space time, there is a band of spinodally 
or parametrically unstable wave vectors, within this band the amplitude of the quantum fluctuations grows. Because 
of the cosmological expansion wave vectors are redshifted into the unstable band and when the wavelength of the 
unstable modes becomes larger than the Hubble radius these modes become classical with a large amplitude and a 
frozen phase. These long wavelength modes assemble into a classical coherent and homogeneous condensate, which 
obeys the equations of motion of the classical inflaton [stUss) . This phenomenon of classicalization and the formation 
of a homogeneous condensate takes place during the first 5—10 e-folds after the beginning of the inflationary stage. 
The full quantum theory treatment in refs. 37, 38] show that this rapid redshift and classicalization justifies the use 
of an homogeneous classical inflaton leading to the following robust conclusions [37i l3^ : 

• The quantum fluctuations of the inflaton are of two different kinds: (a) Large amplitude quantum fluctuations 
generated at the begining of inflation through spinodal or parametric resonance depending on the inflationary 
scenario chosen. They have comoving wavenumbers in the range of e^^~^° lO^^GeV < A: < e^T-&o j^qIS QgV 
and they become superhorizon a few efolds after the begining of inflation. The phase of these long-wavelength 
fluctuations freeze out and their amplitude grows thereby effectively forming a homogeneous classical condensate. 
The study of more general initial quantum states featuring highly excited distribution of quanta lead to similar 
conclusions^^: during the first few e-folds of evolution the rapid redshift results in a classicalization of long- 
wavelength fiuctuations and the emergence of a homogeneous coherent condensate that obeys the classical 
equations of motion in terms of the infiaton potential, (b) Cosmological scales relevant for the observations 
today between 1 Mpc and the Hubble radius had first crossed (exited) the Hubble radius about ^ 50 e-folds 
before the end of inflation within a rather narrow window of about 8 e- folds !l|. These correspond to small 
fluctuations of high comoving wavevectors in the rangeQ g^r-60 iQ^^ GeT^ < k < e^T-60 iq'^o Q^y ^^^qj-q 
Nt > 60 is the total number of efolds. 

• During the rest of the inflationary stage the dynamics is described by this classical homogeneous condensate 
that obeys the classical equations of motion with the inflaton potential. Thus inflation even if triggered by an 
initial quantum state or density matrix of the quantum field, is effectively described in terms of an homogeneous 
scalar condensate. 

The body of results emerging from these studies provide a justification for the description of inflationary dynamics 
in terms of classical homogeneous scalar field. The conclusion is that after a few initial e-folds during which the 
unstable wavevectors (a) are redshifted well beyond the Hubble radius, all that remains for the ensuing dynamics is 
an homogeneous condensate, plus small fluctuations corresponding to modes (b). 

C. Slow roll inflation 

Amongst the wide variety of inflationary scenarios, slow roll inflation|3^^] provides a simple and generic descrip- 
tion of inflation consistent with the WMAP dataflg. In this scenario, inflation is driven by the dynamics of the 
classical coherent and homogeneous condensate of the inflaton field ipit), which obeys the classical equation of motion 

<fi + 3Hip + V'{(p) =0 . (30) 

.2 

Its energy density is given by p{t) = ^ + ^("ys). The basic premise of slow roll inflation is that the potential is fairly 
flat during the inflationary stage. This flatness not only leads to a slowly varying inflaton and Hubble parameter, 
hence ensuring a sufficient number of e-folds, but also provides an explanation for the gaussianity of the fluctuations 
as well as for the (almost) scale invariance of their power spectrum. A flat potential precludes large non-linearities 
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in the dynamics of the fluctuations of the inflaton, which is therefore determined by a gaussian free field theory. 
Furthermore, because the potential is flat the inflaton is almost massless (compared with the scale of V^), and modes 
cross the horizon with an amplitude proportional to the Hubble parameter. This fact combined with a slowly varying 
Hubble parameter yields an almost scale invariant primordial power spectrum. Departures from scale invariance and 
gaussianity are determined by the departures from flatness of the potential, namely by derivatives of the potential 
with respect to the inflaton field. These derivatives are small and can be combined into a hierarchy of dimensionless 
slow roll parameters that allow an assessment of the corrections to the basic predictions of gaussianity and scale 
invariance 18J. The slow roll expansion introduces a hierarchy of small dimensionless quantities that are determined 
by the derivatives of the potential [39ll4ll |: 



The slow roll approximation [sl 139 . l41| corresponds to ~ ?7v ^ 1 with the hierarchy ~ 0{ey), namely ev and rjy 
are first order in slow roll, second order in slow roll, etc. The slow roll variable ey ^ 1 implies that the evolution of 
the inflaton ip is slow, and to leading order (neglecting the second derivatives), the equation of motion 13U|) becomes 

3H{t) LP + V'{if) = (32) 

with 

The second slow roll variable ryy <C 1 implies that the inflationary potential is nearly flat during the inflationary 
stage. During slow roll inflation the number of e-folds, from the time t till the end of inflation, at which the value of 
the inflaton is <pe, is given by 



NMt)] = ^T7Tl 7777^ dLf. (34) 



Mil Kit) y'iv) 

Small fluctuations of the scalar (matter) fields around the classical inflaton lead to small fluctuations in the space-time 
geometry through Einstein's equations. There are two types fluctuations that are relevant: curvature perturbations 
and gravitational waves, both are produced during inflation. The fluctuations of the scalar field generate directly 
fiuctuations in the curvature of space-time, while the expansion of the Universe, itself determined by the dynamics of 
the scalar field generates gravitational waves (see refs.pl Is^). In the linearized approximation different wavevectors 
of the fluctuations evolve differently, the power spectra per logarithmic wave vector interval for curvature (i?) and 
gravitational wave (h) fluctuations arc respectively given by 

/:^l{k)^^{\Ru?) , A2,(fc) = ^(|/i+fep + |/,,,,|2), (35) 

where the Rk are the quanta of curvature fluctuations, /i+,x are the two independent polarizations of the quanta of 
gravitational waves, and the expectation value is in the vacuum state during inflation. These power spectra feature 
power laws in the slow roll regime, 

A|,,(fc) = A|,,(fco) (^- j , A|(fc) = A2(A:o)(^-j , (36) 

namely during slow roll the power spectra of curvature and gravitational wave perturbations are nearly scale invariant. 
The amplitude of the curvature perturbations for wavelengths that re-entered the Hubble radius at the last scattering 
surface are directly related|3,l3 to the temperature anisotropics measured by COBE and WMAPpTl [T5 | -[l8 | . 



A 4 AT 



= (4.67 ±0.27) X 10"^ (37) 

WMAP 



The amplitude and the power laws are determined in slow roll by 

2A.^MUev ' A|(M=^^^^ ' - 1 = "6 + 2 W ; nr = ~2 ey . (38) 
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The WMAP values are for feg = 0.002/Mpc, which yields an upper bound for the scale of the inflationary potential 
during slow roll inflation V^^"^ < 3.3 x lO^^GeV, suggesting a connection between the scale of grand unification and that 
of inflation. While inflation may not be related to phase transitions in the early Universe, the nearly scale invariant 
spectrum of gaussian fluctuations suggests a connection with a critical theory. Ref. '45| provided an effective field 
theory description of inflation akin to the Landau- Ginzburg description of critical phenomena. Slow roll dynamics can 
be organized elegantly in a systematic expansion in powers of l/Ne within a Landau- Ginzburg effective field theory. 
This is achieved by introducing a dimensionless inflaton field and rescaled potential as 



and 



X 



PI 



(39) 



(40) 



where the value of [eq.JSZIl] and eq.ljSHl) fix the scale of inflation M to be M 10^'' GeV and ^ 50 is the 
number of e-folds. To emphasize that the slow roll approximation implies a slow time evolution it is also convenient 
to introduce a stretched (slow) dimensionless time variable r and a rescaled dimensionless Hubble parameter h as 
follows 



t 



H = 



the Einstein-Friedman equation now reads 



^ ' 3 



2 N, 



dx 
dT 



and the evolution equation for the inflaton field x is given by 

1 



d^X , dx 



(41) 



(42) 



(43) 



The slow-roll approximation follows by neglecting the l/N^ terms in eas. (l42ll and (|43|l . Both w{x) and ft.(r) are 
of order N9 for large N^. Both equations make manifest the slow roll expansion as a systematic expansion in 
1/A^e0j Slow roll dynamics, the inflationary scale and all of the observational phenomenology is reproduced 
with X ~ f^(x) ~ C'(l) during inflation. 

As a pedagogical example for slow roll inflation we display the main features of the inflationary dynamics for 
w{x) = (x^ ~ 4)^/16 in figd These figures reveal clearly that inflation ends at r '--^ 10 when /i(r) begins a rapid 
decrease when the field is no longer slowly coasting along near the maximum of the potential x ^ but rapidly 
approaching its equilibrium minimum. The left panel In a(T) depicts the number of e-folds as a function of time [by 
definition Ne{t) = Ina(t)], wavelengths of cosmological relevance today have crossed the Hubble radius during the 
last ~ 50 e-folds before the end of inflation within a narrow window AN^ ~ 8, corresponding to a very small interval 
At ^ I. A successful inflationary scenario requires at least iVg 60 to solve the horizon and entropy problems of 
the standard hot Big Bang 1] which occur in an interval At ~ 3 prior to the end of inflation and during which the 
dimensionless field changes about Ax ^ 1. These are general aspects of a wide range of inflationary scenarii|4^. 
We consider here a translationally and rotationally invariant cosmology where the only source of inhomogeneities 
are (small) quantum fluctuations. Indeed, inhomogeneities cannot be excluded at the beginning of inflation but the 
redshift of scales during inflation by at least ^ 10^^ effectively erases all eventual initial inhomogeneities. 

Inflation is now an established part of cosmology with several important aspects, such as the super horizon origin of 
density perturbations, having been spectacularly validated by WMAP Simple but phenomenologically accurate 
descriptions of inflation invoke an effective field theory for a hornogeneous scalar field, akin to the Landau-Ginzburg 
description of critical phenomena described in section lTlI Al above lial . Analysis of the WMAP data puts very strongly 
pressure on the simple monomial w{x) = G4 X^ at ths Scr level |l8l liof and forthcoming observations of the CMB 
as the Planck satellite and others bear the promise of yielding precise information on the inflaton effective potential 
corresponding to the stage of inflation during which wavelengths of cosmological relevance today first crossed the 
Hubble radius during inflation (see 0,|^Di,|42i| for discussions of the lessons learned so far). 
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V. THE ELECTROWEAK SCALE: PHASE TRANSITIONS AND BARYOGENESIS 



There is a large body of observational evidence that suggests that the there is more matter than antimatter in the 
Universe up to scales of the order of the Hubble radius [43. l45|. The origin of this baryon asymmetry is one of the 
deep mysteries in particle physics and cosmology. The value of this asymmetry is quantified by the ratio 

^ = (44) 

where Ub (ni) is the baryon (antibaryon) density and n-y is the photon density. This is the only free input parameter 
that enters in nucleosynthesis calculations of the primordial abundance of light elements ^3 UM- The agreement 
between the WMAP results[ll El El El and the most recent analysis of the primordial deuterium abundance |46) 
yields 

r;= (6.1±0.3) X 10~^° (45) 

What is the origin of this ratio? namely what is the microscopic mechanism responsible for baryogenesisl . The 
necessary conditions for successful baryogenesis were first identified and outlined by Sakharovl4^: 

• Baryon number violation. 

• C and CP violation: unless these symmetries are violated in any process in which baryons are created and 
annihilated, the rates for baryon production equals that of the reverse reaction and no net baryon can be 
generated by these processes. 

• Departure from equilibrium: in equilibrium the density matrix only depends on the Hamiltonian (and simul- 
taneously commuting operators) which in all microscopic theories is invariant under CPT. Since the baryon 
number operator is odd under CPT, the expectation value of the baryon operator must vanish if the density 
matrix is that of equilibrium^. 



Notably, the standard model of particle physics has the main ingredients for baryogenesis |45l |4 

9 B + L violation in the standard model: Since the weak interactions only involve the left handed quark 
(baryon) and lepton currents there is a quantum mechanical anomaly in their conservation laws|49j: 



Mr 



" 32 7r2 



Wl.W'^" + 5'" B^.B^'^ (46) 



where N f is the number of generations, VF, B are the field strength tensors for the SU{2) and U{1) gauge fields 
and H/^" = e^^"^"^ Wa/s and similarly for B. As a consequence of this anomaly the change in the baryon (and 
lepton) number is related to the change in the topological charge of the gauge field 



AB = Nf ANcs ; Ncs = J d'x e'^'= Tr 



(47) 



W, djWk + ^ ^» Wk 



where Wi is the SU{2) gauge field and Ncs is an integer that characterizes the topological structure of the 
gauge field configuration. Kuzmin, Rubakov and Shaposhnikovlsoj noticed that in the high temperature medium 
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that prevailed in the early Universe there are non-perturbative field configurations, called sphalerons that induce 
transitions between gauge field configurations with different values of Nqs- The sphalerons lead to B+L violating 
processes with transition rate per unit volume estimated to be 51] Tgph ~ ct^ T'^ ln(l/Q;). This estimate for the 
transition rate suggests that sphaleron processes are in thermal equilibrium for 100 GeV < T < 10^^ GeV. 

• CP violation in the standard model Because only left handed quarks and leptons couple to the charged and 
neutral vector bosons that mediate the weak interactions, the standard model violates P maximally. However, 
CP violation is much more subtle and is the result of CP violating phases in the complex Cabibbo Kobayashi 
Maskawa mass matrix for quarks resulting from complex Yukawa couplings to the Higgs field. For Nf generations 
of quarks and leptons there are {Nf — l){Nf — 2)/2 independent phases in the CKM mass matrix, and a non-zero 
value for any of these phases implies CP violation. For Nf > 3 there is at least one CP violating phase, hence 
the standard model, with Nf — 3 does indeed have the possibility of CP violation. Experimentally CP violation 
is observed in the Kq Kq and B B systems. 

• Non-equilibrium: as discussed above weak interaction processes are in LTE down to T ~ 1 MeV, therefore 
the only possibility for non-equilibrium is through a phase transition. Since the expansion rate of the Universe 
is much smaller than the weak interaction rate, it is very likely that a second order phase transition at the 
electroweak scale T ~ 100 GeV would occur in LTE, hence departure from equilibrium requires a strong first 
order phase transition. An estimate of the possibility and strength of a first order electroweak phase transition 
in the standard model is gleaned from a one loop calculation of the effective potential in the SU{2)+ Higgs 
model (i.e, neglecting the [/(I) gauge group) [isf 

V« T) . (f + I + 1=1) (T^ - r.V^ - 1^ + ^ (48, 

where (p = ^J4>*(f>, vq is the Higgs vacuum expectation value at T = and mt is the top quark mass. The 
second term proportional to ip'^ arises from the gauge field contribution and is responsible for a first order 
phase transition. Including the U{1) gauge group changes the above only quantitatively. The one loop effective 
potential as a function of ip has a typical shape as in the right panel in fig. |21 and features a global and a local 
minimum for(45| 

^^^^ = "^"[— + ^- 256^ + ^) 

where niH is the Higgs mass. A measure of the strength the phase transition is the ratio Aip{Tc)/Tc 52] where 
Aip{Tc) is the jump in the order parameter between the two (degenerate) minima at Tc- Successful baryogenesis 
requires this ratio to be > l|52|]. One important aspect that emerges from this simple analysis is that the higher 
the Higgs mass the weaker the phase transition. A study of higher orders [s^ reveals that perturbation theory 
does not give reliable information about the electroweak phase transition for Higgs masses beyond ^ 70 GeV. 
Lattice studies have shown 54] that the ratio A(p{Tc) /Tc < 1 for mn > 45 GeV. // the standard model features 
a strong first order phase transition this phase transition occurs via the formation of nucleated bubbles just 
as described in section IIII Bl Finally, a mechanism for baryogenesis involves transport of baryons through the 
bubble walls as the nucleated bubbles grow and percolate filling the space with the globally stable phase ^Ssj]: 
CP violating interactions of quarks and leptons in the thermal medium with the bubble walls leads to an excess 
of left handed quarks, which sphaleron transitions convert into a net baryon asymmetry. For a discussion of 
these mechanisms see|43 and references therein. 

Caveats: While the standard model features the main ingredients for successful baryogenesis, a substantial body of 
work has revealed that for a Hig gs m ass larger than about 72 GeV there is no first order phase transition but a smooth 
crossover in the standard model|5^ Is^ . The current LEP bound for the standard model Higgs mass mn ^115 GeV, 
all but rules out the possibility of a strong first order phase transition and suggests a smooth crossover from the 
broken symmetry into the symmetric phase in the standard model. While it has become clear that the LEP bound on 
the Higgs mass precludes baryogenesis in the standard model, some supersymmetric extensions of the standard model 
with a stop lighter than the top may be able to explain the observed baryon asvmmetrv |56l IstI Iss f . An alternative 
scenario for baryogenesis proposes that a primordial asymmetry between leptons and antileptons or leptogenesis is 
responsible for generating the baryon asymmetry ^48. 59] . The leptogenesis proposal depends on the details of the 
origin of neutrino masses and remains a subject of ongoing study. 
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VI. THE QCD PHASE TRANSITION IN THE EARLY UNIVERSE 

One of the most spectacular epochs m the early Universe is the QCD transition, when quarks and gluons become 
confined in hadrons. In the early Universe the baryon asymmetry is very small [see eq. (|45ll ]. and at RHIC (and soon 
at LHC) it is expected that the mid-rapidity region for central collisions is also baryon free (see the discussion in 
section IVIII Al below') . For the purpuse of cosmology, one would wish to directly measure Tqqd and the equation of 
state from experiments with relativistic heavy ions, but it turns out not to be a simple task. The problem arises from 
extremely different time scales: 10~^ s in the cosmological QCD transition (very close to thermal equilibrium), but 
only 10^^^ s in the laboratory (out-of-equilibrium effects may be important). 

Since QCD is asymptotically free, it is expected that at high temperature a perturbative evaluation of the equation 
of state in terms of a weakly interacting gas of quark and gluons should be reliable. However near the hadronization 
phase transition the nature of the degrees of freedom changes from quarks and gluons to hadrons and QCD becomes 
non-perturbative. The only known first principle method to study QCD non-perturbatively in a wide temperature 
range is lattice gauge theory (LOT). Over the last several years there has been steady progress in the study of the 
QCD phase diagram with and without chemical potential including light and heavy quarks p30ll6lL l6^. As explained 
above, the problems to incorporate a finite chemical potential in LGT (see [g^l) is of no concern to cosmology and the 
mid-rapidity region of central collisions at RHIC and LHC, because the relevant baryochemical potentials /is/T ^ 1. 

For an extended review of the cosmological consequences of the QCD transition and the physics of the first second 
of the Universe see [b^ . 

A. The QCD transition and equation of state: 

1. Lattice gauge theory results 

It has been established that lattice QCD without dynamical quarks (quenched approximation) exhibits a thermal 
first-order phase transition 0| at a critical temperature of T^, w 270 MeV ^E^. This is also in agreement with the 
expectation of a first-order phase transition from the simplest bag model (SSl EE (see below) . Unfortunately, 
the situation is not clear for dynamical quarks and is especially unclear for the physical values of the up, down and 
strange quark masses. 

For dynamical quarks, lattice QCD calculations provide a range of estimates for Tc. In the case of two-flavour QCD 
Tf, sa 175 MeV 68), 6^, whereas for three-flavour QCD ~ 155 MeV [b^, almost independent of the quark mass. 
For the most interesting case of two light quark flavors (up and down) and the more massive strange quark, a value of 
Tc ~ 170 MeV has been obtained recently, both from standard 70] and improved [tJ] staggered quarks. Accordingly, 
in the discussion that follows we adopt a transition temperature Tc = 170 MeV, keeping in mind that the systematic 
uncertainty is probably of the order 10 MeV. 

The order of the phase transition and the value of Tc is still under investigation. For massless quarks the theoretical 
expectation is a second order transition for two quark flavors and a first-order transition for three and more quark 
flavors [72] . On the lattice, for two ligh t quarks the results are inconclusive. The predicted universality class is not 
confirmed so far (see discussions in [60|); most recent studies even claim to find hints for a first-order transition [tJ. 
For three flavours close to the chiral limit the lattice results clearly indicate that the phase transition is of first order 
(see 1131), as expected from theory. Some older simulations suggested that this holds true for the physical case as well 
The latter result was obtained using the Wilson quark action, whereas results with standard |70j and improved 
Tlf staggered quarks indicate a crossover for the physical quark masses. 

A microscopic description of phase transitions in QCD requires also a reliable assessment of the equation of state 
(EoS). Furthermore, as it will be described below, a hydrodynamic description of the space-time evolution in relativistic 
heavy ion collisions also requires knowledge of the equation of state at high temperature. The consensus that seems 
to be emerging is that for the physical masses of two light (up and down) and one heavier (strange) quark there is a 
sharp crossover between a high temperature gas of quark and gluon quasiparticles and a low temperature hadronic 
phase without any thermodynamic discontinuities. This is displayed in fig. |5] which summarize results from LGT for 
the energy density and pressure (both divided by to compare to a free gas of massless quarks and or gluons) as a 
function of T/Tc |62l|. The most recent simulations of the equation of state are reviewed in [g^. 

These figures clearl y re veal a sharp decrease in the energy density and pressure at T = Tc, the value of the energy 
density is £c ~ 6 T^ |62j | thus predicting an energy density Ec 0.7GeV/fm^. Furthermore the high temperature 
behavior is not quite given by the Stephan-Boltzmann law, see fig. (jSJ, suggesting that even at large temperatures 
the plasma is not described by free quarks and gluons up to temperatures T ^ 4.0 Tc ~ 700 Mev. 



22 




In the hydrodynamic limit, the speed of sound 

(50) 



de 



s 



is a quantity of central interest. A strong decrease in the speed of sound, already above Tc, has been observed in 
lattice QCD, see 75, 76] for the most recent data. Consistent with previous results from quenched QCD and two-flavor 
QCD, c^(Tc) « 0.1 when approaching the critical temperature from above (in the QGP). The right panel of fig. |H1 
displays the speed of sound Wd\ clearly showing a dramatic decrease for T < 2Tc and approaching 1/3 for T Tc 
in agreement with an ultrarelativistic gas of quarks and gluons. 



2. The bag EOS and a first order phase transition 

If the transition from T > Tc to T < Tc is continuous but sharp as evidenced by the lattice data in fig. |H| the 
behavior may not be too different from an actual transition which may be modelled by a simpler (EoS) which would 
allow an analytic treatment and to be included in a hydrodynamic description. The bag model 65] provides a semi- 
phenomenological description of an (EoS) that features a quark-hadron transition. The simplest version considers the 
thermodynamics in two different regions: for T > Tc a gas of massless quarks and gluons and confinement is accounted 
for in terms of a bag constant Bg5], for T < Tc a, gas of free massless pions (in the chiral limit). At T = Tc quarks, 
gluons and pions coexist in equilibrium at constant pressure and temperature Pc,Tc. In this model the pressure and 
energy density for vanishing chemical potential are given by 

p> = 5> ^ - B ; ey^gy^T* + B ; for T > , (51) 

P<^g<—T^ ; £< = g<—T'^ ; forT<Tc, (52) 

where B is the bag constant with a typical value B^ ^ 200 MeVjs^. For a description of a quark-hadron transition 
g> = 37, 47.5 for Nf = 2,3 respectively and g< — 3 for a gas massless pions. The two regions are joined together by 
a flat coexistence curve with constant pc, Tc according to a Maxwell construction (see sec lIII B 3|) . The value of Tc is 
determined by this coexistence condition and is given by 

. ^ ,53) 

For Nf — 2 and three massless pions Tc ~ 145 Mev which is not too far from the lattice results Tc ~ 170 MeV. 
The flat coexistence region describes a mixed phase in which a quark-gluon plasma is in coexistence with a hadron 
gas with a vanishing isothermal sound speed, just like in the liquid-gas transition. The proportion of each phase in 
coexistence is determined by the lever-rule described in section IIII B 31 The bag equation of state is therefore given 
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by 

p^^-^B ; for e>e>{T,) , p = i ; for s < e^i^) , (54) 

with a flat coexistence hne for e<(Tc) < p < e^{Tc). This (EoS) is displayed in fig. O compare to the liquid-gas 
equation of state in the left panel of fig. 

Since the Helmholtz free energy density is f{T) = e — t s = —p with s the entropy density, it is clear that for T > Tc 
the free energy of the quark gluon phase is smaller and for T < Tc the hadron phase has a smaller free energy. Since 
on the coexistence line p — Pc, T = Tc are constant, there is a jump in the entropy, namely a latent heat between the 
QGP and the hadronic phase, hence the bag model (EoS) leads to a first order phase transition with a latent heat 

l = Tcl\s = e>{Tc)-e<{Tc)=AB , (55) 

independent of the number of degrees of freedom. In the quenched approximation on the lattice, the latent heat was 
determined to be Z « 1.4 T^ [zSlzll; which is a factor of 5 smaller than expected from a quenched bag model. 
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FIG. 9: Bag (EoS) for a massless pion gas for T < Tc and quarks and gluons in a bag for T > Tc- 



More realistic versions of the bag equation of state have been constructed in ref. |66', . The importance of the 
bag (EoS) not only stems from its simplicity and the fact that it models a sharp crossover as evidenced by the lattice 
data but also because it is widely used in the hydrodynamic approach to ultrarelativistic heavy ion collisions [see the 
discussion in section IVIII A II below]. 

Besides the latent heat, the surface tension is the crucial parameter for the nucleation of bubbles in a first-order 
phase transition. The surface tension <j = {dW /dA)y is the work dW that has to be done per area dA to change the 
phase interface at fixed volume. The absence of surface excitations in hadronic spectra suggests that cr^/^ <C B^^^ 
[79] in the bag model. A self-consistent calculation of the surface tension within the bag model shows that the surface 
tension vanishes for massless quarks and gluons if no interactions besides the bag constant (i.e. as = 0) are taken into 
account [l^- This can be cured by introducing short-range interactions |^ or by including the strange quark mass 
[S^,!^- For the surface tension rather small values are found from quenched lattice QCD, a « 0.015 [tsLIs^ . 



B. Cosmological consequences of the QCD transition 



In the radiation dominated Universe for T > Tc the particle content is: quarks, gluons, leptons and photons. It is a 
good approximation to treat all particles with to <C 3 T as massless. Above the QCD transition ^quarks = (7/8)12A''/ 
(number of quark flavours) and ggiuons = 8, below the QCD transition ghadrons quickly drops to zero, i.e. pions disappear 
at r - 40 MeV. Therefore at the QCD epoch g> = 51.25 (61.75) , 5< = 17.25(21.25) without (with) strange quarks 
and without (with) kaons, respectively. At the QCD transition the Hubble radius is about 10 km, corresponding to 
scales of 1 pc or 3 light-years today. The mass inside a Hubble volume is (47r/3) e{Tc) R% - IA/q. This mass IS 
redshifted oc (1 -I- z) as the Universe expands, because it is made up of radiation. The mass of cold dark matter in 
a comoving volume is invariant, M^dm = [(1 + ^cqj/i^ + z)] M{z). At the QCD transition (1 + Zcq)/(1 -I- z) ^ 10"^ 
(recall z^q ^ 10'*), and thus 

M§°^ - 10~^Mq . (56) 

The Hubble time at the QCD transition, tgcD ~ 10~^ s, is much longer than the relaxation time scale of particle 
interactions. Thus, leptons, photons, the QGP and the hadron gas (HG) are in thermal and chemical equilibrium 
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during cosmological time scales. All components are in LTE and form a radiation fluid. There conserved quantum 
numbers for this radiation fluid, lepton and baryon numbers (in the standard model). However, the corresponding 
chemical potentials are much smaller than the temperature at the QCD epoch. In this situation we can use the 
equation of state calculated in thermal lattice QCD and apply it directly to study the equilibrium aspects of cosmology. 
In standard cosmology, all lepton chemical potentials are assumed to vanish exactly. The baryon number density 
rtB follows from the ratio of baryons to photons at Big Bang Nucleosynthesis (BBN), rj = (nB/?^7)BBN, and the 
conservation of baryon number and entropy, ne/s constant: 

nB(Te) = 77 s{Tc) = (8.2 ± 0.4) x 10"" s(Te) , (57) 

V S / BBN 

using the value of rj given in ea. l|45|) . Finally, the baryon number inside a Hubble volume is given by 
at the beginning of the transition and about twice that value at the end^. 



1. A first-order QCD transition 

If the cosmological QCD transition is of first order, it proceeds via bubble nucleation "ss", "sl*!. Its typical duration 
is O.ltQCD- From the small values of surface tension and latent heat found in lattice QCD calculations [s^l, the 
amount of supercooling is found to be small fsJl. Hadronic bubbles nucleate during a short period of supercooling, 
A^sc ^ 10"'^ tqcD- In a homogeneous Universe without 'dirt' the bubbles nucleate owing to thermal fluctuations 
(homogeneous nucleation), with a typical bubble nucleation distance of ^g^J 

dnuc - 1 cm ~ 10"^ Rn . (59) 

The hadronic bubbles grow very fast, within lO^^tqcD, until the released latent heat has reheated the Universe to 
Tc- By that time, just a small fraction of volume has gone through the transition. For the remaining 99% of the 
transition, the HQ and the QGP coexist at the pressure Phg(^c) = Pqgp(^c)- During this time the hadronic bubbles 
grow slowly and the released latent heat keeps the temperature constant until the transition is completed. The energy 
density decreases continuously from eQGp(T'c) at the beginning of the transition to ehgCTc) when the transition is 
completed. 

In first-order phase transitions that we know from our everyday experience, like the condensation of water drops 
in clouds, the drops (bubbles) are nucleated at impurities ('dirt'). This could happen in the early Universe as well. 
Candidates for cosmic 'dirt' are primordial black holes, monopoles, strings, and other kinds of defects. Of course, 
the existence of any of these objects has not been verified so far. Under these circumstances the typical nucleation 
distance may differ significantly from the scenario of homogeneous nucleation. For the QCD transition, heterogeneous 
nucleation has been studied in some detail in ref. 86]. 

When the magnitude of primordial temperature fluctuations is of the same order or larger than the typical super- 
cooling, the transition proceeds via inhomogeneous bubble nucleation. The mean nucleation distance results from the 
scale and amplitude of the temperature fluctuations. 

Homogeneous nucleation : The probability to nucleate a bubble with critical radius (i.e. the minimal bubble size 
that can grow after its formation) by a thermal fluctuation per unit volume and unit time is given by (see discussion 
in Sec. HI B.2) 

I{T) = Io{T) cxp(^-^) , (60) 

with AWc = 16 TT cr^/[3 (phg — Pqgp)^]- For dimensional reasons Iq ^ C T^, with C = 0{1). A more detailed calcula- 
tion of lo within the bag model has been provided in ■ It was shown in Ref. that the temperature dependence 
of the prefactor Jo can be neglected for the calculation of the supercooling temperature Tgc in the cosmological QCD 
transition. 



^ This formula is correct if no black holes are formed during the QCD transition and if the quark nuggets that might have formed evaporate 
before the BBN epoch. 
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FIG. 10: Sketch of a first-order QCD transition via homogeneous bubble nucleation: above the critical temperature the 
Universe is filled with a quark-gluon plasma (Q). After a small amount of supercooling the first hadronic bubbles (H) nucleate 
at ti, with mean separation cinuc- At t2 > ti these bubbles have grown and have released enough latent heat to quench the 
formation of new bubbles. The supercooling, bubble nucleation, and quenching takes just 1% of the full transition time. In the 
remaining 99% of the transition time the bubbles grow following the adiabatic expansion of the Universe. At the transition 
is almost finished. The shrinking QGP drops are separated by the typical distance dnuc- 



For small supercooling A = 1 — T/Tc ^ 1 we may evaluate (phg ~ Pqgp){T) by using the second law of thermo- 
dynamics, i.e. phg — Pqgp ~ ^ a, and thus 

/(A)«/o(T,) exp(-^/A2), (61) 

with A = 16 n cr^/(3 P T^) and /o(rc) ~ T^- Note that this result does not depend on the details of the QCD 
equation of state. For the values of / = 1.4 and cr = 0.015 from quenched lattice QCD ^ A ^ i y. 10"^. In 
the bag model ^ « 5 x 10"^ {a/T^ f. 

The amount of supercooling that is necessary to complete the transition. Age, can be estimated from the schematic 
case of one single bubble nucleated per Hubble volume per Hubble time, which is 



0(A„ 



A 



-,1/2 



4 \n{T,lHQCD) 



4 X 10"^ (62) 



for the values of I and a from quenched lattice QCD. For the bag model we assume a < 0.1 T^, which implies that 
Asc < 6 X lO-'^. 

The time lapse during the supercooling period follows from the conservation of entropy and reads 

Atsc/tQCD = A,c/(3 cl) = 0(10-3) . (63) 

Here we used the relation = dlns/dlnT for the speed of sound in the supercooled phase. For realistic models 
< Cs{A) < 1/V3, and in the bag model Cs(A) = l/VS. 

After the first bubbles have been nucleated, they grow most probably by weak deflagration [s^. Issl Issl 189^. The 
deflagration front (the bubble wall) moves with the velocity Wdefl ^ l/v^ 0|. The energy that is released from the 
bubbles is distributed into the surrounding QGP by a supersonic shock wave and by neutrino radiation. This reheats 
the QGP to Tc and prohibits further bubble formation. Since the amplitude of the shock is very small [s^, on scales 
smaller than the neutrino mean free path, heat transport by neutrinos is the most efficient. Neutrinos have a mean 
free path of 10~^i?H at T^. When they do most of the heat transport, heat goes with Wheat — C>{c). For larger scales, 
heat transport is much slower. Figure ^| shows a sketch of the homogeneous bubble nucleation scenario. 

Let us now calculate the mean bubble separation, dnuc and the final supercooling. Age, for a scenario with weak 
deflagration. Bubbles present at a given time have typically been nucleated during the preceding time interval 

Atnuc = I/{dI/dt) . (64) 

Using the relation between time and supercooling, dA/dt = 3 c^/tqcD, we find 

Atnnc/tQCD = Al/{6Acl)^O{10-^) and A^uc = A,e = 0(10-^) A,c . (65) 
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During the time interval Atnuc each bubble releases latent heat, which is distributed over a typical distance « 
2 viiea.t A^nuc- This distance has a weak dependence on the precise value of Age, but the bubble nucleation rate 
increases strongly with A until one bubble per volume (A^nuc Whcat)^ is nucleated. Therefore the mean bubble 
separation is 

dnuc « 2 Wheat A^c « ^ ^ -^H = ©(IQ-*' Rh) - O(lcm), (66) 

where we used Wheat — 0(0.1), 3 = 0(0.1), which gives a typical value for the nucleation distance. The suppression 
of bubble nucleation due to already existing bubbles is neglected. 

The estimate ea. (|66|) of the mean bubble separation applies if the released latent heat by means of sound waves and 
by neutrino free streaming is sufficient to reheat the QGP to T^, i.e. to quench the nucleation of new bubbles. On the 
other hand the typical bubble separation could be given by the rate of release of latent heat, i.e. by the bubble wall 
velocity Wdcfl- Since the period of supercooling lasts about 1% of the time needed for completing the entire first-order 
phase transition, 1% of the QGP must be converted to HQ in the process of sudden reheating to Tc, the bubble radius 
at quenching must therefore reach 0.2 of the bubble separation, -Rbubbic ~ 0.2 dnuc- With i?bubbie ~ VdcR ^tauc, and 
using the above relation d^uc ~ 2 Uhoat A^nuc, we require Wdcfl > 0.4?;heat for consistency. If Wdcfl is smaller than this, 
the limiting factor for quenching is the rate of release of latent heat by bubble growth, and the bubble separation is 



2 Wdefl Atnuc ~ ^ ^ , (67) 

6 ci A 



i.e. the bubble separation will be smaller than the estimate in eq. (|66|l . 

We are now in a position to improve the estimate of Asc- one bubble nucleates in the volume (wheat Atnuc)'^ during 
A^nuc- This can be written as 

1 « (Wheat Atnuc)^ Atnuc/(tsc) , (68) 

which in terms of the supercooling parameter Age is given by: 
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Alt exp --nr N 10^^ A^^,^ exp — — . (69) 



While the pre-exponential factor is smaller than the naive estimate ea. H62f) by a factor of 10^°, the amount of super- 
cooling is just 20% larger than in ca. l|^ . i.e. Age = 5 x 10^^, confirming, that that numerical prefactors in ea. l|^ 
are irrelevant in the calculation of Age. 

In summary, the scales on which non-equilibrium phenomena occur are given by the mean bubble separation, which 
is about 10^^i?H- The entropy production is tiny, i.e. AS/S ^ 10^^, since the supercooling is small ^ 10"'^. After 
supercooling, which lasts 10"'^ tgcD, the Universe reheats in Atnuc ~ 10~^ tqcD- After reheating, the thermodynamic 
variables follow their equilibrium values and bubbles grow only because of the expansion of the Universe. 

Inhomogeneous nucleation: The local temperature T{t, x) of the radiation fluid fluctuates, because cosmological 
perturbations have been generated during cosmological inflation's^. Let us denote the temperature fluctuation by 
At = 5T/T. Inflation predicts a Gaussian distribution of perturbations (see sectiorUvjl : 



P{At) dAr = —p^ exp 
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If one allows for a tilt in the power spectrum of density fluctuations, the rms temperature fluctuation reads [see ea. H36|) 
and ref.|i3], 

A5.-^ « 10-4 (3 ^2)3/4 ^ (71) 

where fco is the wave number of the mode that crosses the Hubble radius today. For the scale-invariant spectrum 
{us = 1) it is foundlij] A™''(A:ood) ~ 2 x 10"^ 

The detailed analysis in ref.[9ll leads to the conclusion that the picture of homogeneous bubble nucleation, where 
bubbles form from statistical fluctuations, is false for the most probable cosmological scenarios. 

A new scenario for the cosmological QCD transition was developed in ref.joj in which a small scale cut-off in the 
spectrum of primordial temperature fluctuations from coUisional damping by neutrinos [o^ . l93j | is included. This 
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FIG. 11: Sketch of a first-order QCD transition in tlie inhomogeneous Universe (from [oH'): at ti the coldest spots (dark 
grey) are cold enough to render the nucleation of hadronic bubbles (H) possible, while most of the Universe remains in the 
quark-gluon phase (Q). At t2 > t\ the bubbles from the cold spots have merged and have grown to bubbles as large as the 
fluctuation scale. Only the hot spots (light grey) are still in the QGP phase. At ts the transition is almost finished. The last 
QGP drops are found in the hottest spots of the Universe. The mean separation of these hot spots can be much larger than 
the homogeneous bubble nucleation separation. 

study revealed that at the QCD transition neutrinos travel freely on scales A,y_mfp ~ 10^® i?H, and that fluctuations 
on the diffusion scale of neutrinos are washed out by the time of the QCD transition: 

Ai/-diff = 2 V A^-mfp c tqcD ~ 10 i?H • (72) 

Thus the old picture of homogeneous bubble nucleation still applies within the small homogeneous patches of Agmooth — 
IQ-^ i?H- 

The compression time scale for a homogeneous patch is 5t — Agmooth/cs ^ 10^^ t-QCD- If the compression time 
scale is larger than Ainuc the temperature fluctuations are frozen with respect to the time scale of nucleations. 

A sketch of inhomogeneous bubble nucleation is shown in fig. 1111 The basic idea is that temperature inhomogeneities 
determine the location of bubble nucleation. In cold regions, bubbles nucleate first. In general we have two possible 
situations: 

1. If Anuc > A™*', the temperature inhomogeneities are negligible and the phase transition proceeds via homoge- 
neous nucleation (see section IVI B ip . 

2. If Aiiuc < A™*^, the nucleation rate is inhomogeneous and we have to consider the scenario sketched in fig. 1111 

A first attempt to analyse inhomogeneous nucleation has been given in l9l|. According to 91], the nucleation distance 
dnuc exceeds the scale Agmooth, if 

Asmooth < 2 Jp^/\}™^ ^jj, (73) 

3 cj 

If A™ > 5 X 10~^, it is quite likely that this condition is met. In that case we can conclude that the typical 
inhomogeneity scale in the baryon distribution is inherited from the scale of density inhomogeneities in the radiation 
fluid at the end of the QCD transition. The effect in terms of length scales is at least two orders of magnitude larger 
than the nucleation distance in homogeneous nucleation and is 0(1 m), which is of interest for inhomogeneous BBN. 

2. Effects from a first-order QCD transition 

Let us now briefly summarize the effects that have been suggested to emerge from the cosmological QCD transition. 
There are two kinds of effects: the effects that have been found in the mid 80s and early 90s stem from the bubble 
scale and they thus affect scales A < dnuc- The formation of quark nuggets, the generation of isothermal baryon 
fluctuations, the generation of magnetic flelds and gravitational waves belong to the effects from the bubble scale. 
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In recent years it was found that there is another class of possible consequences from the QCD transition, which 
are connected to the Hubble scale and therefore affect scales A < i?H- Among these effects are the amplification of 
inhomogeneities and later formation of cold dark matter clumps, the modification of primordial gravitational waves, 
and the enhanced probability of black hole formation during the QCD transition. 

Quark nuggets/Strangelets: In 1971 Bodmer joj] suggested the possibility that strange quark matter might 
be the ground state of bulk matter, instead of ^^Fe. Later Witten [2^ rediscovered this idea. Strange quark matter 
was further studied by Farhi and Jaffe [t^- The idea of strange quark matter is based on the observation that the 
Pauli principle allows more quarks to be packed into a fixed volume in phase space if three instead of two flavours are 
available. Thus the energy per baryon would be lower in strange quark matter than in nuclei. However, the strange 
quark is heavy compared with up and down quarks, and this mass counteracts the advantage from the Pauli principle. 
No strange quark matter has been found experimentally so far [OQj . The issue of stability of strange quark matter has 
not been settled yet; for a recent review see [o^ . 

Witten [q^ pointed out that a separation of phases during the coexistence of the hadronic and the quark phase 
could gather a large number of baryons in strange quark nuggets These quark nuggets could contribute to 

the dark matter today or affect BBN [T^ ll,'^ |. At the end of the transition the baryon number in the quark 

OOP 

droplets could exceed the baryon number in the hadron phase by several orders of magnitude, could be close 

to nuclear density. However, it was realized that the quark nuggets, while cooling, lose baryons. The quark nuggets 
evaporate as long as the temperature is above ~ 50 MeV |99j . Quark nuggets may survive this evaporation if they 
contain much more than ~ 10^*^ baryons initially Llj3d |. This number should be compared with the number of baryons 
inside a Hubble volume at the QCD transition, which is 10''*. Thus, the mean bubble nucleation distance should be 
> 3 X 10"^i?H ^ 300 m so as to collect enough baryons. This seems impossible from todays perspective, as explained 
above. 

Inhomogeneous nucleosynthesis: Applegate and Hogan (10 1| found that a strong first-order QCD phase transi- 
tion induces inhomogeneous nucleosynthesis. It is extremely important to understand the initial conditions for BBN, 
because many of our ideas about the early Universe rely on the validity of the standard (ho moge neous) BBN sce- 
nario. This is in good agreement with observations |l2Lll3| . In inhomogeneous nucleosynthesis |l02| . large isothermal 
fluctuations of the baryon number (the remnants of the quark droplets at the end of the QCD transition) could lead 
to different yields of light elements. As a minimal requirement for an inhomogeneous scenario of nucleosynthesis , the 
mean bubble nucleation distance has to be larger than the proton diffusion length, which corresponds to ^ 3 m 103| 
at the QCD transition. This is two orders of magnitude above recent estimates of the typical nucleation distance 8q|. 

Although values for i] dramatically different from those in the standard BBN are excluded both from measurements 
of the light element abundances and from the CMB, it might be possible to alter the primordial abundance of heavy 
elements {A > 7) in inhomogeneous scenarios |l04j . 

Cold dark matter clumps: Scales A that are of the order of the Hubble radius i?H are not sensitive to details of 
the bubbles. It was reported in Refs. |93 . ll05| | that the evolution of cosmological density perturbations (see Sec. HI) 
is strongly affected by a first-order QCD transition for subhorizon scales, A < Ru- 
in the radiation-dominated Universe subhorizon density perturbations perform acoustic oscillations. The restoring 

1/2 

force is provided by pressure gradients. These, and therefore the speed of sound Cs = {dp/de)g (on scales much 
larger than the bubble separation scale) drop to zero at a first-order QCD transition |105| . because both phases coexist 
at the pressure Pc only (a is the scale factor of the Universe): 



de(a)/da 

It stays zero during the entire transition and suddenly rises back to the radiation value Cs = l/\/3 after the transition. 
A significant decrease in the effective speed of sound Cg during the cosmological QCD transition was also pointed out 
by Jedamzik ^IQQj. 

As the speed of sound drops to zero, the restoring force for acoustic oscillations vanishes and density perturbations 
for subhorizon modes fall freely. The fluid velocity stays constant during this free fall. Perturbations of shorter 
wavelengths have higher velocities at the beginning of the transition, and thus grow proportional to the wave number 
k during the phase transition. The primordial Harrison-Zel'dovich spectrum of density perturbations is amplified on 
subhorizon scales. The spectrum of density perturbations on superhorizon scales, A > i?H, is unaffected. At T ~ 1 
MeV the neutrinos decouple from the radiation fluid. During this decoupling the large peaks in the radiation spectrum 
are wiped out by coUisional damping '92]. 

Today a major component of the Universe is dark matter, most likely CDM. If CDM is kinetically decoupled from 
the radiation fluid at the QCD transition, the density perturbations in CDM do not suffer from the neutrino damping. 
This is the case for primordial black holes or axions, but not for supersymmetric dark matter. At the time of the 
QCD transition the energy density of CDM is small, i.e. ecdm(?c) ~ 10^* erad(?c)- CDM falls into the potential wells 
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provided by the dominant radiation fluid. Thus, the CDM spectrum is amphfied on subhorizon scales. The peaks in 
the CDM spectrum go non-linear shortly after radiation-matter equality. This leads to the formation of C DM clumps 
with mass < 1O^^"M0. Especially the clumping of axions has important implications for axion searc hes |l07j . If the 
QCD transition is strong enough, these clumps could be detected by gravitational femtolensing |108| |. 



3. Damping of gravitational waves at the QCD transition 

In principle, primordial gravitational waves (e.g. from cosmological inflation) present a cle an pr obe of the dynamics 
of the early Universe, since they know only about the Hubble expansion. As was shown in 'l09l| a step is imprinted 
in the spectrum of primordial gravitational waves by the cosmological QCD transition. This step does not allow us 
to tell the difference between a first-order transition and a crossover, but an estimate of the transition temperature 
and a measurement of the drop in effective number of relativistic degrees of freedom would be possible. 

Primordial gravitational waves are predicted to be generated during inflation |ll(]| | and could be detected by observ- 
ing the so-called B-mode (parity odd patterns) polarisation of the CMB. Inflation predicts an almost scale-invariant 
energy density per logarithmic frequency interval for the most interesting frequencies (~ 10~^ Hz for pulsar timing, 
~ 10^"^ Hz for LISA, and ~ 100 Hz for LIGO and VIRGO) of the gravitational waves. The energy fraction in 
gravitational waves, per logarithmic interval in frequency /, is defined by 



f^gw(/) = / 
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Figure IT^ shows the transfer function rigw(/)/r2gw(/ ^ fqco) from the cosmological QCD transition. The typical 
frequency scale is 
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which corresponds to the mode that crosses the Hubble horizon at the end of the bag model QCD transition. Scales 
that cross into the horizon after the transition (l.h.s. of the figure) are unaffected, whereas modes that cross the 
horizon before the transition are damped by an additional factor « 0.7. The modification of the differential spectrum 
has been calculated for a first-order (bag model) and a crossover QCD transition. In both cases the step extends over 
one decade in frequency. The detailed form of the step is almost independent from the order of the transition. 
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FIG. 12: The modificatio n of t he energy density, per logarithmic frequency interval, for primordial gravitational waves from 
the QCD transition (from |109|I 1. 



The size of the step can be calculated analytically [109j. Comparing the differential energy spectrum for modes 
that cross into the horizon before and after the transition gives the ratio 



^gw(/ > fqcp) 
^gwif < Iqcd) 



0.696 



(77) 



for the QCD transition. 

In fig. 1121 we i ndic ated the frequency range (^ ^ , 
timing residuals |lll| . Unfortunately, today's technology does not enable us to detect primordial gravitational waves 
at frequencies around 10~^ Hz, because their expected amplitude is too small. 



1 yr ^) in which limits on VL„^ have been reported from pulsar 
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VII. OBSERVING THE EARLY UNIVERSE: THE EARLIEST IS THE LATEST. 

As discussed above, during the inflationary stage the acceleration of the scale factor is positive, while in the 
radiation or matter dominated eras it is negative. Ea. H27|) shows that wavelengths grow faster than the Hubble radius 
during inflation but slower during radiation and matter dominated eras. Perturbations of cosmological relevance with 
wavelengths between our Hubble radius today to galactic scales exited the Hubble radius during a window of ~ 5 — 8 e- 
folds about 55 e-folds before the end of inflation thus becoming causally disconnected from microphysics. Wavelengths 
relevant to CMB anisotropies re-entered at the time of recombination prior to photon decoupling. Therefore CMB 
anisotropies provide direct information about the inflationary stage. Contrary to this case, the wavelengths of quantum 
fluctuations produced inside the Hubble radius during the radiation or matter dominated eras are always inside the 
Hubble radius and are causally influenced by microphysical processes. Hence observable consequences of processes 
that occurred during these eras are generally indirect. The baryon asymmetry is perhaps a remnant of the the 
microphysics of the standard model (or extensions thereof) and perhaps the QCD phase transition left an observable 
footprint as discussed in sees I VI B2l and IVI B 31 The time scale between the QCD phase transition ~ lO^secs and 
that of nucleosynthesis ^ 200 sees is very large on the scale of strong interactions ~ 10~^^ sees which tend to erase 
any potentially interesting signature. Thus CMB observations, 400000 years after the Big Bang, provide a window 
into the earliest phenomena that took place during inflation, while microphysical processes within the standard model 
of particle physics occur much later, are causal all throughout the evolution, and interactions blur their observable 
signals. An important possible remnant of these phase transitions are primordial magnetic fields. 

A variety of astrophysical observations including Zeeman splitting, synchrotron emission, Faraday rotation mea- 
surements (RM) combined with pulsar dispersion measurements (DM) and polarization measurements suggest the 
presen ce of larg e scale magnetic fields 112. 113j . The strength of typical galactic magnetic fields is of the order 
~ /i ^ 11121 ITTsf an d the y are correlated on very large scales up to galactic or even larger reaching to scales of cluster 
of galaxies ~ 1 MpcfTT^. The origin of these large scale magnetic fields is still a subject of discussion. It is currently 
agreed that a variety of dynamo mechanisms ar e effi cient in amplifying seed magnetic fields with typical growth 
rates F ^ Gyr~^ over time scales ~ 10 — 12 Gyr |112| . There are different proposals for the origin of the initia l seed 
at different stages in the history of the early Universe, in particular during inflation, and or ph ase t ransitions jll2j . 
Primordial (hyper) magnetic fields may have importa nt co nsequences in electroweak baryogenesis jH^ , Big Bang nu- 
cleosynthesis (see 112|), the p olarization of the CMB [iT5I | via the same physical processes as Faraday rotation, and 
structure formation ^ll2j . |llfil |. If the electroweak and/or t he ch iral phase transitions occurred out of equilibrium they 
could be a significant source of primordial magnetic fields llF^ . Thus cosmic magnetic fields may be one of the few 
observational relics of primordial phase transitions. The study of the origin of cosmic magnetic fields is one of the key 
science projects of the forthcoming Square Kilometer Array (SKA) 118]. 

VIII. STUDYING PHASE TRANSITIONS WITH ACCELERATORS 
A. Ultrarelativistic heavy ion collisions: seeking the Quark Gluon Plasma 

The program of ultra relativistic heavy ion collisions (URHIC) whose primary goal is to study the phase diagram 
of QCD began almost two decades ago with the fixed target heavy ion programs at the AGS at Brookhaven and the 
SPS at CERN. A summary of the results of these efforts mainly through the Ph -\- Ph experiments at SPS-CERN 
provided compelling evidence in favor of the existence of a new state of matter |2^. The program continues at the 
Relativistic Heavy Ion Collider (RHIC) at BNL. Unlike the fixed target experiment SPS at CERN, RHIC is a collider 
experiment which currently studies Au + Au collisions with center of mass energy -y/s ~ 200AGev and luminosity 
~ lO^^cm^^s^^. The future ALICE (A Large Ion Collider Experiment) heavy ion program at LHC is expected to 
study Pb + Pb collisions with cm. energies up to \/s ~ 5ATev and luminosities ~ lO^^cm^^s^^. In these collisions 
the heavy nuclei can be pictured in the CM frame as two Lorentz contracted pancakes. For Au + Au collisions the 
size of each pancake in the direction transverse to the beam axis is about 7 fm. At RHIC and LHC energies most of 
the baryons are expected to be carried away by the receding pancakes (the fragmentation region) while in the region 
of the collision a large energy (density) is deposited in the form of quark pairs and gluons. At least two important 
mechanisms for energy deposition in the collision region are at work IHtI lllfll Il2nt Il2lj : i) the establishment of a 
strong color electric field (flux tube) that eventually breaks up into quark-antiquark pairs when the energy in the field 
is larger than the pair produc tion t hreshold and ii) the partons (quarks and gluons) inside the colliding nuclei interact 
and redistribute their energy |122| . An estimate of the energy deposited in the collision region has been provided by 
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Bjorken [122 



1 dEx 



(78) 



with To ~ Ifm/c, Ra ~ 7fm for Au and dEx/dy is the transverse energy per unit r apidi t y wh ich is measured. Parton- 
parton scattering is expected to lead to thermahzation on time scales ^ Ifm/c |ll9l Il20| . Af ter t he quark-gluon 
plasma achieves LTE, the evolution is conjectured to be described by hydrodynamic expansion |124| . As the QGP 
expa nds a nd cools, the temperature falls near the critical temperature and the confinement and chiral phase transitions 
occur [T25II . Upon further cooling the quark-gluon plasma hadronizes, the hadrons rescatter until the hadron gas is 
dilute enough that the mean free path is larger than the mean distance between hadrons. At this point hadrons freeze- 
out and stream out freely to the detectors from this last scattering or freeze- out surface. This picture is summarized 
in fig. El below. Estimates based on this picture and on detailed numerical evolution (126ii] suggest that at RHIC the 
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FIG. 13; Heavy Ion Collisions. 



quark gluon plasma lifetime is of order ~ 10/m/c while the total evolution until freeze-out is ~ 50 — 100/m/c. 



1. Hydrodynamics, lattice gauge theory and the equation of state. 



Although the evolution of the quark gluon plasma from the initial state described by the parton distribution of the 
colliding nuclei until freeze-out of hadrons clearly requir es a non- equilibrium description, a hydrodynamic picture of 
the evolution is both useful and experimentally relevant |124L Il27j . In the approximation in which quarks and gluons 
are strongly coupled in the sense that their mean free paths are much smaller than the typical wavelength for the 
variation of collective phenomena, the QGP can be described as a fluid in LTE, for which the energy momentum 
tensor is of the form 

T^'^ = {e+p)u^ u'' -p g'"' , (79) 

with e, p the energy density and pressure respectively, g'^'^ is the metric and 

= 7(x) (l,i/(a;)) (80) 

is a local 4-velocity vector with 7(0;) the local Lorentz contraction factor. 

In this description the dynam ical evolution of the fluid is obtained from the conservation laws of energy- momentum, 
baryon number and entropy |6^ll2lLll24lll25| | with an equation of state p = p{e) to close the set of equations. A simple 
and phenomenologically useful description is provided by Bjorken's model for longitudinal expansion 67, 121, 12J5, 12JI- 
There are three important ingredients in Bjorken's models: i) the central rapidity region where the energy deposited 
and the formation of the QGP takes place well separated from the fragmentation region which is the region of 
the receding pancakes where most of the baryons are (see fig. I13|l . ii) The thermodynamic variables are invariant 
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under boosts along the beam (longitudinal) axis, this is based on the observation that the particle distributions are 
invariant under these boosts in the central rapidity region and iii) the expansion only occurs along the beam axis, 
i.e, longitudinal expansion. In the baryon free region, local thermodynamics implies that the entropy density s in the 
local rest frame of the fluid is related to the energy density e and pressure p as 



e + p ^ T s . 

It is convenient to introduce the space time rapidity t] and proper-time r as 



^ = iln 



t + z 



t - z 



(81) 
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with z the coordinate along the beam axis (longitudinal) and the fluid rapidity 9 by writing the local velocity eq. H8U|) 
as — (cosh 6', 0, 0, sinh0). Bjorken's model assumes the fluid to be composed of free streaming particles for which 
Vz = z/t in which case the fluid rapidity 9 becomes the space-time rapidity n . Boost invariance along the longitudinal 
direction entails that e,p, s, T are all functions of proper time onlv jl23tll24j . 

The energy-momentum conservation equation d^jT^'^ leads to two equations by projecting along the direction 
and perpendicular to it using the projector g^" — u^u'^. Under the assumption that the central region is ba ryon free, 
the conse rvation of entropy and energy and momentum lead to the following equations (for details see Il2li Il23i 
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Assuming an equation of state p{t) = p{e{T)) and combining eas. H83() with the thermodynamic relation H81() (for 
baryon free plasmas) one finds that the evolution of the temperature is given by 



r dT 2 2 dP 

Which for constant speed of sound results in the cooling law 
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To 



(85) 



Obviously, the form of eq. H83|l is similar to the energy conservation equation JHJ in cosmology with the expansion rate 
a/a = l/(3r) when the proper time r is identified with the comoving time t in the cosmological setting. The similarity 
becomes even more remarkable for the case of the QGP being modelled as a radiation fluid (which is expected to be a 
good approximation at high temperature) since in this case = 1/3 and the connection with a radiation dominated 
cosmology with scale factor a{t) cx is evident. 

To find the general evolution equations for the QGP, an equation of state is needed. It is at this stage where 
the connection with lattice gauge theory (LOT) is made. LGT obtains the thermodynamic functions in equilibrium 
and these are input in the hydrodynamic descr iptio n as local functions of space-time under the assumption of LTE. 
However, practically hydrodynamic simulations [123 ^^e a bag (EoS). For a discussion of the parameters in the bag 
(EoS) as well as a justification for using this equation of state which features a strong first order phase transition 
instead of the lattice (EoS) which features a crossover see [l27j . 



B. Predictions and observations pre-RHIC 

Several experimental signatures had been associated with the formation of the QGP |ll9l ll28| , and the heavy ion 
program at SPS focused on several of them. We summarize the observables that have been proposed as telltales of a 
QGP and the data gath ered by the SPS from Pb + Pb and Pb + Au collisions, which taken together provide a hint 
of evidence [H HI [llS [ill for a QGP, although many of them could have alternative explanations. 



1. J /'^ suppression and strangeness enhancement 



The J/^f is a narrow cc bound state, which, if produced in the early stages of the collision can probe the QGP 
because its lifetime is longer than that of the QGP and decays into dilepton pairs which leave the plasma without 
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scattering. The original suggestion [iSOj is that when the screening length of the color force is smaller than the size 
of the cc bound state, this narrow resonance will melt. This argument suggests that charmonium suppression could 
provide evidence for a QGP. Alternative lv a suppression could result from scattering with hard gluons in the plasma 
and the dissociation of the boimd state 130]. A normal suppression of charmonium is expected on the grounds that 
once formed the cc bound state interact with other nuclcons inside the nucleus. This expected suppression is studied 
in proton nuclcon collisions and extrapolated to nucleon-nuclcon collisions. This is considered the normal suppression 
in contrast to the abnormal suppression expected from the presence of a plasma. The left panel of fig. 1141 shows the 
data gathered by the NA50 collaboration at the SPS-CERN 131]. 

This figure reveals an abnormal suppression when the energy density is e > 2.5GeV/fm^. The en ergy density in 
this figure has been computed with Bjorken's formula ea. H78|) . The NA50 collaboration at CERN-SPS |l32j combined 
data for J/^* suppression from the NA38 and NA50 experiments. The analysis reveals that while for the most 
peripheral (largest impact parameter) collisions the suppression can be accounted for by nuclear absorption, there is 
no saturation in the suppression in the most central Pb + Pb collisions and that the observed suppression pattern can 
be understood in a deconfinement scenario. This report concluded that the J/vP suppression pattern observed in the 
NA50 data provides significant evidence for deconfinement of quarks and gluons in Pb + Pb collisions. Strangeness 
enhancem ent a long with chemical equilibration are some of the earliest proposals for clear signatures of the formation 
of a QGP |l33l | . The main idea is based on the estimate that the strangeness equilibration time in a hot QGP is of the 
same order as the expected lifetime of the QGP (~ 10/m/c) produced in nucleus-nucleus collisions. Two important 
aspects of this estimate make strangeness enhancement a prime candidate: if strangeness attains chemical equilibrium 
in the QGP, this equilibrium value is significantly higher than the strangeness production in nucleon-nucleon collisions. 
Also s tran geness production through hadronic rescattering or final state interactions was estimated to be negligibly 
small jl33j| . In the QGP, color deconfinement leads to a large gluon density that leads to the creation of ss pairs, 
furthermore chiral symmetry makes the strange quark lighter thus lowering the production threshold. This situation 
is in contrast to the case of hadronic rescattering or final state interactions where the production of pairs of strange 
quarks has large thresholds and small cross sections . The usual measure of strangeness enhancement is through 
the ratio As = 2{ss)/ {uu + dd). The right panel in fig. E| displays as a function of -/s for nucleon-nucleon as well 
as nucleus-nucleus colHsions {S + S, S + Ag, Pb + Pb) at SPS. The data displayed in fig. ^^is taken as evidence that 
nucleus-nucleus collisions at SPS are creating a hot state of matter. 
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FIG. 14: Left panel: anomalous J/^ suppression as a function of the initial energy density. From |l3ll ]. Rig ht pa nel: The ratio 
\s as a function of ^/s, the energy of the collision for nucleon-nucleon and nucleus-nucleus collisions. From |l34jl 



2. Electromagnetic probes: dileptons and direct photons 

Elect romagnet ic probes: e'^'e" or fJ.~^ fJ-~ dilepton pairs and direct (prompt) photons are prime probes of the hot 
plasma jTH ll^. since once produced they leave the plasma without further interactions because their mean free path 
is much larger than the typical size of the plasma. The p vector meson is important because it decays into dileptons 
and its lifetime is ^ 1/to/c. Therefore once it is produced in the hadron gas it decays within the hot hadronic plasma 
and the produced dileptons carry direct information from the plasma. Thus while dileptons produced from the decay 
of the p meson do not yield evidence of the earlier stages in the QGP, they do offer information on the hadronic 
stage. The left panel in fig. E| presents the data gathered by the CERES- NA45 collaboration for the invariant mass 
spectrum for electron-positron pairs from 158 AGeV Pb+Au collisions at the SPS-CERN. The solid line represents the 
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expected spectrum from the decays of hadrons produced in proton-nucleon and proton-proton collisions extrapolated 
to Pb + Au collisions and is the sum of the contributions shown in the graph. There are two remarkable features 
in this graph: a clear enhancement of dileptons in the region 2 MeV < Mg+g- < 700 MeV and that instead of the 
p meson peak at mp = 770 MeV there is a broad distribution. The excess of dileptons in the small invariant mass 
region cannot be explained by charged pion annihilation .23.] . It is remarkable that the dilcpton enhancement is below 
the putative p peak and that there is no hint of the p at 770 MeV!. The current understandin g o f these fea tures is 
that the medium effects result in a shift in the p meson mass as well as a change in its width |2j, 112(1 11351 . Thus 
while this interpretation does not directly yield information on the QGP, it does support the picture of a hot gas of 
hadrons, mainly pions which is the main interaction channel of the p vector meson. Direct photons are conceptually 
a clean direct probe of the early stages of the QGP. Photons are produced in the QGP by several processes: gluon- 
to-photon Compton scattering off (anti)quark q{q)g — * 9(9)7 and quark-antiquark ann ihilation to photon and gluon 
qq g"f and to the same order (see Kapusta et. al. and Aurenche et. al. in '1365) (anti)quark bremsstrahlung 
1l{9) ~^ 1Q{9)i and quark-antiquark annihilation with scattering qqq{g) q(,g)l- Detailed calculations including 
screening corrections |l36j reveal that direct photons from the QGP could provide a signal that could be discriminated 
against the signal from the hadronic background. This work indicates the theoretical feasibility of direct photons as 
direct probes of the early stages of the QGP. The WA98 collaboration at SPS-CERN reported their analysis for the 
first observation of direct photons from Pb + Pb collisions with ^/i = 158 AGeV |l37j |. Their data is summarized in 
the right panel of fig. 1151 The transverse momentum distribution of direct photons is determined on a statistical 
basis and compared to the background photon yield predicted from a calculation of the radiative decays of hadrons. 
The most interesting result is that a significant excess of direct photons beyond that expected from proton-induced 
reaction at the same y/s is observed in the range of transverse momentum greater than about 1.5 GeV/c in central 
collisions. A com paris on of the data to the theoretical predictions for direct photons from an equilibrated QGP was 
performed in ref |l38j . The conclusions are that while it is not clear if SPS reached the energy density to form the 
QGP, the data indicates a hot and dense phase that could be its precursor. It has recently been shown |l39j | that 
nan- equilibrium effects in an expanding QGP formed in RHIC collisions lead to an enhancement of direct photons in 
the region of transverse momentum pt ^ 2 GeV. 




FIG. 15: Left figure: invariant mass spectrum of e^e" pairs from Pb + Au collisions at 158 AGev at the SPS-CERN 
(CERES/NA45). From 140]. Right figure: ratio of total measured yield of direct photons to hadronic background vs. 
transverse momentum from Pb + Pb at 158 AGev for peripheral and central collisions. From |l37j| . 



3. Quantum spinodal decomposition and pion domains 



For massless up and down quarks, QCD has an SU{2)l ® SU{2)fi symmetry which is spontaneously broken down 
to SU{2)l+r, this is the chiral phase transition. It has been argued that for massless quarks the equilibrium aspects 
of this transition is described by the universality class of the 0(4) Heisenberg ferromagnet_30j. The Landau-Ginzburg 
effective field theory for this universality class is given by the effective potential 

y($)-a(T-rg)$2 + g($2)2 . $^(^^^)^ (86) 

where a ipj^ip is the pseudoscalar chiral order parameter field that acquires an expectation value below the 
critical temperature, and tt ~ t/jj^rip (t" = isospin matrices) is the triplet of pions, which are the Goldstone bosons 
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associated with the symmetry breaking. In an URHIC, chiral symmetry is expected to be restored if a QGP is 
formed with T > T^. Upon expansion and rapid coohng chiral symmetry is again broken. // the chiral phase 
transition occurs out of equilibrium, the ensuing quench can bring the hadronic phase into the spinodal region of the 
effective potential (18611 . As a result, quantum spinodal instabilities would lead to the formation of correlated domains 
inside which the order parameter is disoriented resulting in a disorient ed ch iral con(ie7isate(DCC),3^ |141||. Spinodal 
instabilities would result in the formation and growth oi pion doma,ins \l4'^ with distinct observational signatures in 
the pion distribution. A systematic study by the WA98 collaboration at SPS-CERN'143'] measured charged particle 
multiplicities in 158 AGeV Pb + Pb collisions and ruled out the formation of (DCC) domains at 90% CL. The Minimax 
collaboration studied the possibility of (DCC)-like events at Fermilab's Tevatron in pp collisions 144] and found no 
statistically significant signals of large isospin fluctuations. 



C. News from RHIC: 



The four experiments at RHIC have prese nted their analysis of the first three years of operation at the top energy 
-y/s = 200 AGeV for Au + Au collisions jl45j . The combined bod y of results not only confirms most of the findings 
at SPS but also reveal novel and unexpected phenomenafl^, ■ To begin with, the results on particle multiplicity 
when combined with Bjorken's estimate for the value of the energy density eq.jTHJ yield e ~ 5.5GeV/fm^ for tq ~ 
Ifm/c. This energy density is almost one order of magnitude larger than the critical value obtained from (LGT), 
Elgt ~ 0.7GeV/fm^ (see section (|VI A therefore these experiments probe QCD in the deconfined regime. At least 
two important results stand out: 



1. Elliptic flow, ideal hydrodynamics and early thermalization 



If a QGP is formed early in the collision the pressure drives a collective expansion or flow. An important measure of 
this collective flow is elliptic flowQA^ in collisions with non- vanishing impact parameter (non-central). If the reaction 
plane is taken to be the x — z plane, with z the beam axis, the triple differential momentum distribution for hadron 
h is expanded in a Fourier series in terms of the angle 4> in the x — y plane as 

'^^'^ '^^^ [l + 2wi(y,pT)cos</) + 2«2(y,PT)cos2(/)+---] . (87) 



dydp^d(j) 2iT dydp'j^ 

The first coefficient vi vanishes at midrapiditv |l27|| . the second W2, the ellip tic flow parameter is a measure of the 
hydrodynamic flow in peripheral collisions driven by pressure gradients| l27| . A remarkable result from RHIC ex- 
periments is that ideal hydrodynamic calcu lations of V2 are in excellent agreement with the combined experimental 
results from STAR [14^ and PHENIX |148| | for a many hadron species, up to ~ 2 GeV as shown in fig. ^1 A 




FIG. 16: Combined STAR and PHENIX data on V2. FromlH. 



hydrodynamic calculation requires two inputs: an initialization (proper) time tq at which LTE is assumed to hold 
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and an equation of state. The initial stage that leads to LTE requires knowledge beyond the realm of applicability 
of hydrodynamics. Hydrodynamic calculations just input initial data at a given (proper) time hypersurface. In most 
hydrodynamic approaches a bag-like (EoS) is used[l27] with typical values for the bag constant B^^^ ~ 200 — 250 MeV 
which yield a critical temperature Tc ^ 160MeV. Such equation of state yields a strong first order transition, with 
a latent heat I ~ IGeV/f m'^. While such strong first order phase transition is in contradiction with the (EoS) from 
(LGT), it is argued in ref. |l27| that the latent heat is consistent with the (gradual) change in entropy density across 
the transition in the lattice data. Agreement between the hydrodynamical calculations and the experimental values 
of V2 requires a very short initialization time scale tq 0.6fm/c (127.] and systematic studies of different equations 
of state in ref.ji^ points out that the data is reproduced if the (EoS) includes a soft point at which the speed of 
sound becomes anomalously small as featured by a quark-hadron mixed phase. As t he p lasma enters the mixed phase 
pressure gradients become very small and anisotropies in the transverse flow stall |l27j . Furthermore an important 
aspect of the hydrodynamic results is that ideal hydrodynamics yields a remarkable agreement with the experimental 
data on collective flow. Taken together these results suggest the following: 

• Early thermalization: the data is reproduced if the initialization proper time at which LTE is assumed to hold 
is To ~ 0.6fm/c. 

• Ideal hydrodynamics: the data is reproduced by ideal hydrodynamics, namely a strongly coupled fluid with a 
short mean free path and small viscosity coefhcients. 

• Phase transition: the study of different (EoS) in ref.[6^ conclude that hydrodynamics with an (EoS) that feature 
a phase transition is consistent with the elliptic flow data, while an (EoS) without a phase transition is not. 

The emerging consensus is that the new state of matter formed in URHIC is a strongly correlated state or strongly 
coupled quark gluon plasma, a strikingly different situation from the original conjecture of a QGP of weakly interacting 
quark and gluon quasiparticles. Ideal hydrodynamics combined with such early thermalization scale leading to a 
rapid equilibration can only result from strong correlations. However, and alternative interpretation has been recently 
argued: that the RHIC data, in particular the elliptic flow parameter V2 can be interpreted from an incomplete 
equilibration [l49( | . Clearly the understanding of the new findings at RHIC is still evolving. 

2. A new initial state? 

A new picture of the in itial state prior to the onset of LTE and hydrodynamical evolution is emerging |2j|. Current 
theoretical understanding ^15^ suggests that the ground state of high energy nuclei prior to the collision is akin to a 
non-perturbative Bose-Einstein condensate of gluons, a color glass condensate. Evidence for a rapid growth in the 
gluon distribution function as a function of x, the fraction of the lo ngitudinal momentum carried by the gluon was 
first revealed by deep inelastic scattering experiments at HERA p^^ . Gluon self-interaction results in a characteristic 
momentum scale below which the gluon distribution function saturates, at RHIC this scale is Qs ^ IGeV. The 
density of gluons below this scale is dNg/d\n{l/x) ~ l/a^ which is very large because the strong interaction coupHng 
as becomes small at this scale. The evidence for this novel initial state at RHIC emerges from the suppression or 
quenching of jets and has been observed in comparing d + Au and Au + Au collisions 22, 150J. This high gluon 
density initial state is conjectured to be the precursor of a QGP^| and it is possible jl52j that this highly correlated 
and non-perturbative initial state is responsible for the early thermalization required by the hydrodynamic approach. 
These novel aspects of URHIC will continue to be the focus of the experimental and theoretical programs. 

D. Little Bang vs. Big Bang 

One of the stated goals of the URHIC program was to open a window to the early Universe by creating the 
conditions that prevailed when it was about 10 /isecs old. While there are similarities between URHIC and early 
Universe cosmology, there also are differences which are even more pronounced by the recent results from RHIC. 

Similarities: i) The early Universe is almost baryon free, and the entropy is dominated by radiation, in central 
collisions at RHIC (and the forthcoming LHC) the central rapidity region is nearly baryon free, ii) Ideal hydrodynamics 
and longitudinal expansion are very similar to the Hubble expansion of the cosmic fluid. 

Differences: i) In the early Universe at the time of the QCD phase transition, the expansion time scale 1/H ^ 
10~^ sees was much longer than the time scale of strong interactions ^ 10~^^ sees, hence the crossover behavior 
revealed by the lattice data implies that the QCD phase transition(s) (deconfinement-confinement and chiral symmetry 
breaking) occurred in LTE. In URHIC, the expansion time scale is of order 10~^^secs, ten times the time scale of strong 
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interactions, and non-equilibrium effects are more likely, ii) Since in the early Universe cooling occurs via Hubble 
expansion from an initial temperature or energy density many orders of magnitude larger than the energy density 
required for the QCD phase transition ^ IGeV/fm"^, undoubtedly a weakly interacting gas of quarks and gluons, 
namely a weakly interacting QGP was present as a consequence of asymptotic freedom. Contrary to this expectation, 
the results from RHIC seem to suggest a strongly interacting liquid, a very different state of matter. Furthermore, the 
initial condition in the early Universe, prior to the QCD scale is that of a radiation dominated fluid in LTE slowly 
cooling via adiabatic expansion. RHIC (and HERA) reveal a novel initial state of strongly correlated gluons with a 
non-perturbative 'color glass condensate'. Whether the strongly coupled quark gluon liquid is a consequence of such 
strongly correlated initial state and whether we can learn aspects of QCD during the first few /xsecs after the Big 
Bang will continue to be the focus of the experimental programs at RHIC and LHC. 



E. The nuclear liquid gas phase transition: spinodal decomposition and nuclear fragmentation. 

The interaction between nucleons in nuclei is similar to the Van der Waals potential which acts between molecules: 
a short distance repulsive core and a long distance attractive tail. For this reason is was suggested that the nuclear 
interaction should lead to a liquid-gas phase transition in nuclei [l53l |. This original work suggested that if the equation 
of state of nuclear matter is of the Van der Walls type, collision experiments may bring excited nuclei into the spinodal 
region of the phase diagram in which spinodal instabilities towards phase separation would lead to the spectacular 
fragmentation of nuclei. Thi s obs ervation launched considerable theoretical efforts to yield a better understanding of 
the nuclear equation of state'l54ll| and several experimental facilities to study this phenomenon. The most commonly 
used nuclear effective interaction to study the thermodynamics of nuclear matter is of the S ky rme t-y p el 1 5 5 1 . More 
modern approaches to the nuclear equation of state are in terms of effective nuclear field theories jl56lll57| of hadrons 
and meso ns where the hadronic interaction results from the exchange of tt, cr, p and w mesons. At the mean field 
level |l57j the hadronic interaction does feature a short range repulsion and a long range attractive tail that yield 
an equation of state with a fiat region of coexistence akin to the Maxwell construction in a Van der Walls equation 
of state. These theoretical studies suggest that the critical temperature for nuclear matter is Tc ~ 18 — 20 MeV. 
From the experimental point of view a number of accelerator facilities: GANIL, GSI, LNL, MSU, Dubna have 
been developed to study the physics of nuclear multifragmentation as a consequence of the nuclear liquid gas phase 
transition. These experimental facilities study low energy he avy ion collisions involving several types of nuclei. For 
example heavy ion reactions of Xe + Sn at 32 MeV/nucleon |l58l| are studied at GAN IL, a nd 25 — 35 MeV/nucleon 
Au + Cu colfisions at the K1200-NSCL Cyclotron at Michigan State University (MSU)'l59|. The main experimental 
signature of the nuclear liquid-gas phase transition is expe cted to be the decay of highly excited nuclei through the 
process of multifragmentation via spinodal instabilities [T60|. Expe rimentally the process of multifragmentation leads 
to abnormally large fiuctuations in the distribution of nuclei |l6l| . these fluctuations becoming more pronounced as 
the critical temperature is approach ed. The d istribution of fragments should favor a particular size, associated with 
the most spinodally unstable modes flS^. Il6lj . In finite nuclei the transition is smoothed by the finite size and the 
signatures are not expected to be as sharp as predicted by bulk thermodynamics in an infinite medium. Theoretical 
modelling of the dynamics is based o n transpo rt or quantum molecular dynamics simulations 1 1 54| . There is now a 
substantial body of experimental data'l58', '159] that seems to confirm the presence of such large fiuctuations with an 
impressive identification capability thanks to high performance Att de tectors. A summary of the experimental results 
from the INDRA detector at GANIL has been recently Dresented |l62j| . Some of the remarkable experimental evidence 
reported by the INDRA collaboration are'l62!|: 

• Statistically significant evidence for spinodal fragmentation in Xe+Sn collisions at energies 32—45 MeV/ nucleon. 

• Bimodality: in coexistence the order parameter features bimodality from the two coexisting phases. Such 
bimodality was reported in the distribution of fragments in Au + Au collisions at energies in the range 60 — 
100 MeV/nucleon. 

• Negative specific heat: the microcanonical heat capacity is a thermodynamic quantity that features a negative 
branch in the coexistence region, signaling again an instability. The results of the experiments are in agreement 
with such a negative branch of the microcanonical heat capacity for finite nuclei in a region of coexistence. 

These experimental findings seem to validate the picture of the nuclear liquid-gas phase transition and pave the way 
towards a detailed understanding of the nuclear equation of state in a regime of temperatures and density of relevance 
for understanding nuclear matter in supernovae explosions and the emerging neutron stars. As more high precision 
experiments are carried out at these facilities a more detailed picture of the equation of state of nuclear matter will 
emerge, as well as a first, spectacular confirmation of the dynamics of a phase transition in nuclear matter. 
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IX. SUMMARY, CONCLUSIONS AND OUTLOOK: 



This article presents an excursion of the early and the present Universe exploring the standard models of cosmology 
and particle physics. The focus is on the equilibrium and non-equilibrium physical aspects and potential observational 
consequences of the phase transitions predicted by particle physics. General relativity, statistical mechanics and the 
standard model of particle physics lead to a detailed understanding of the evolution of the Universe and to time- 
marks at which phase transitions may have occurred. A detailed introduction to static and dynamical aspects of 
phase transitions with relevance to cosmology and the experimental program that studies the phase diagram of quark 
and nuclear matter is provided. 

An important aspect of the connection between particle physics and cosmology is the theoretical prediction of an 
early stage of inflation during which the visible Universe expands almost exponentially. This inflationary stage solves 
outstanding problems of the standard Big Bang model and has been spectacularly confirmed by CMB measurements. 
Inflation is now an integral part of the standard model of cosmology and inflationary dynamics has a prominent 
role in the evolution of the Universe. We have discussed the potential observable consequences of the electroweak 
and the QCD phase transitions in the early Universe and highlighted the fact that while anisotropies in the CMB 
provide direct information of the inflationary phase which is the earliest, observable consequences of the electroweak 
and QCD phase transition, which took place much later are indirect. The URHIC program initiated with the AGS 
at BNL, continued with SPS and CERN and currently with RHIC at BNL has provided a wealth of evidence in 
support of a novel state of matter, the QGP. We summarize the experimental evidence from SPS and RHIC and 
discuss some of the most recent discoveries at RHIC. While the results from all four experiments at RHIC provide 
compelling evidence for a new state of matter, the evidence suggests a strongly correlated quark-gluon liquid instead 
of a weakly interacting gas of quarks and gluons which presumably existed in the early Universe prior to the QCD 
phase transition. We have also discussed theoretical and experimental studies of the liquid-gas phase transition in 
cold nuclear matter. Recent results from several experimental programs offer compelling evidence for a liquid-gas 
transition and spinodal instabilities in nuclear matter leading to multifragmentation of excited nuclei. These studies 
yield a detailed knowledge of the equation of state of nuclear matter in a regime of temperature and density relevant 
of supernovae and neutron stars. 
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